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CHAPTER  1 

INTRODUCTION  AND  HISTORICAL  REVIEW 

1.1  Introduction 

When  a  transparent  material  is  mechanically  deformed  its 
index  of  refraction  is  affected.  The  physical  basis  of  the  refrac¬ 
tive  index  change  was  first  considered  by  Pockels,  who  gave  a 
phenomenological  description^  ^ .  In  1922  Brillouin  predicted  that 
an  ultrasonic  acoustic  wave  propagating  in  a  transparent  material 

would  give  rise  to  a  periodic  modulation  of  the  material's  refrac- 

(2) 

tive  index  which  would  diffract  an  incident  light  beam,  and 

(3) 

Debye  and  Sears  observed  the  effect  ten  years  later. 

More  recently,  the  advent  of  the  laser  has  provided 
renewed  interest  in  the  acousto-optic  interaction  of  both  an 
experimental  and  theoretical  nature.  Devices  for  laser  modula¬ 
tion,^^  beam  steering, ^ ^  electronic  filtering, ^  dye  laser 

( 7 )  ( ft  ^ 

tuning,  digital  signal  multiplexing  and  demultiplexing,' 

spectrum  analysis,  ^  and  correlation ^ 10 ’  have  been  demonstrated 

and  analyzed,  to  name  a  few  representative  applications. 

In  the  next  section  we  introduce  some  simple  models 

which  are  useful  in  describing  several  features  of  the  interaction. 

Subsequent  sections  treat  acousto-optical  materials  and  the 

acousto-optic  interaction  in  surface  acoustic  waves. 


2 


1.2  Acousto-Optic  Interactions:  Raman-Nath  and  Bragg  Scattering 

For  a  single  frequency  ultrasonic  wave  traveling  in  the  z- 
direction  in  some  isotropic,  homogeneous  material,  we  can  repre¬ 
sent  the  refractive  index  in  the  region  of  the  wave  as 


n(z,t)  =  n  +  n,  sin  (oj  t  -  k  z) 
o  1  a  a 


(1.1) 


where  nQ  is  the  unpertubed  refractive  index  and  n^  is  the  maximum 
change  in  refractive  index  arising  from  the  presence  of  the  sound 
wave.  In  terms  of  the  elastic  strain  S,  the  elasto-optic  constant 
p,  and  the  unperturbed  refractive  index,  n^  is  given  by^1-^ 


n 


1 


1  : 
2  no 


P  s 


(1.2) 


The  velocity  of  the  acoustic  wave,  v  =  0)  /k  ,  is  normally  about 

Q*  &  & 

five  orders  of  magnitude  less  than  the  velocity  of  light  in  the 
material.  Consequently  an  optical  beam  incident  on  the  sound 
wave  effectively  encounters  a  static  phase  grating  of  periodicity 
given  by  the  acoustic  wavelength  A.  It  is  well  known  that  a 
spatially  and  temporally  coherent,  monochromatic,  collimated  light 
beam  normally  incident  on  such  a  phase  grating  is  diffracted  in 
the  far  field  into  several  orders , ^ ^ ^  shown  schematically  in 
Fig.  1.1.  The  angle  of  the  nth  order  from  the  optic  axis,  0n,  is 
given  by 


n  A 


sin  8  = 

n 


of 


(1.3 


in  which  \  ^  is  the  free-space  wavelength  of  the  incident  light. 
The  intensity  of  the  n^*1  order  diffracted  light,  normalized  to 


k 


that  of  the  incident  light,  is  given  by 


r  -  <♦>• 


(l.M 


Jn  is  the  ntJl  order  Bessel  function  and  <j>  is  the  phase  retardation 
introduced  on  the  light  due  to  the  presence  of  the  sound  wave. 


■  2tt  t 

*  =  T  nl  L 


(1.5) 


in  which  X  is  the  wavelength  in  the  medium  and  L  is  the  inter¬ 
action  length.  The  interaction  length  is  the  distance  over  which  the 
light  and  sound  overlap,  measured  along  the  direction  of  optical 
propagation . 

This  description  of  the  acousto-optic  interaction  in 

which  the  acoustic  wave  is  considered  as  a  phase  grating  was 

(12) 

first  given  by  Raman  and  Nath.  Under  certain  circumstances 

to  be  discussed,  the  effect  of  the  acoustic  wave  is  better 
described  by  considering  the  wave  as  a  series  of  partially 
reflecting  mirrors  corresponding  to  the  maxima  of  the  refractive 
index  wave,  Eq.  (1.1).  In  analogy  to  X-ray  diffraction,  when  the 
light  is  incident  at  the  Bragg  angle  0_.  given  by 

*.  ■  -It  •  m 

the  reflections  constructively  interfere  in  one  direction  giving 
rise  to  a  single  diffraction  order  of  intensity 


X  .  2  ,1 

—  sin  (-  <p) 

o 


(1.7) 


5 


Bragg  diffraction  is  shown  schematically  in  Fig.  1.2(a), 
and  in  Fig.  1.2(b)  is  the  momentum  diagram  for  Bragg  diffraction. 
The  momentum  diagram  corresponds  to  a  particle  model  of  the  inter¬ 
action  in  which  a  photon  and  phonon  are  annihilated  with  the 

(13) 

simultaneous  creation  of  a  new  photon.  Conservation  of 

momentum  provides  the  result  that  the  scattered  light  exits  at 
an  angle 


0 


(1.8) 


and  conservation  of  energy  provides  that  the  frequency  of  the 

scattered  light  is  upshifted  with  respect  to  the  incident  light 

by  the  acoustic  frequency.  The  shift  in  frequency  could  also  have 

been  arrived  at  by  considering  the  doppler  frequency  shift  of 

(lU) 

the  light  wave  reflected  from  moving  mirrors. 

A  unified  theory  treating  both  Raman-Nath  and  Bragg 
scattering  has  been  given  by  Klein  and  CoolJ1'^  in  which  it  is 
shown  that  the  quantity  Q,  defined  by 


Q 


2t r 


Xof  L 


n 


(1.9) 


o 

describes  the  nature  of  the  acousto-optic  interactior.  For 
Q  «  1,  Raman-Nath  diffraction  occurs.  For  Q  »  1,  Bragg  diffrac¬ 
tion  occurs.  In  the  region  Q  'v  1  an  intermediate  type  of  diffrac¬ 
tion  occurs. 

The  results  of  this  section  are  developed  more  precisely 
within  the  framework  of  electromagnetic  theory  in  Chapter  2,  for  the 


Fig.  1.2  Bragg  Diffraction 

a)  Acoustoi-optic  interaction  with  bulk  wave 

b)  Momentum  conservation  for  Bragg  interaction 
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special  case  of  surface  acoustic  waves.  In  the  following  section 
a  discussion  of  acousto-optic  materials  is  given. 


1.3  Acousto-Optic  Material  Considerations 


The  acousto-optic  interaction  can  he  observed  in  both 

liquids  and  solids,  in  the  latter  case  in  amorphous  and  crystalline 

(16) 


media  belonging  to  all  symmetry  classes 
materials  it  is  useful  to  define  certain  figures  of  merit. 
The  most  widely  used  are 


In  comparing  different 
(17,18,19) 


Ml  - 


M2  * 


7  2 
n  '  P 
o 

P  Va 

n6P2 

o 

P  v  3 


(1.10) 


(1.11) 


and 


M_ 


7  2 
n  !  P 
o 

P  v  2 


(1.12) 


in  wich  p  is  the  material’s  mass  density.  The  quantity  M2  is  a 
material  constant  which  is  a  measure  of  diffraction  efficiency. 
In  particular,  the  phase  retardation  <J>  appearing  in  Eqs.  (l.U) 
and  (1.7)  can  be  written  in  the  form 


IT 


(§)M2 


p 

ac 


(1.13) 


in  which  the  ratio  (D/H)  is  the  inverse  aspect  ratio  of  the  acoustic 
beam  waist,  and  P  is  the  acoustic  power.  The  figure  M.  takes  into 

cLC  i. 

account  diffraction  bandwidth  as  well  as  efficiency,  as  does  in 


in  the  special  case  that  the  acoustic  beam  height  is  as  small  as 
the  optical  beam  diameter.  The  relevancy  of  these  figures  is  dis¬ 
cussed  in  detail  by  Maydan.^*^ 

(21) 

Pinnow  has  proposed  guidelines  for  the  selection  of 

acousto-optic  materials  which  provide  a  qualitative  estimate  of 

the  figures  of  merit  of  a  material  from  its  chemical  formula  and 

density.  His  conclusion  is  that  there  is  a  good  possibility  of 

finding  materials  with  high  figures  of  merit  among  the  oxides  which 

are  dense  and  relatively  soft  and  contain  heavier  cations.  This  is 

so  because  softness  is  related  to  low  sound  velocity,  density  with 

high  refractive  index,  and  heavy  cations  with  higher  photoelastic 

constants.  His  guidelines  have  led  to  the  discovery  of  at  least 

(22) 

two  new  acousto-optic  materials,  Te02  and  FbMoO^. 

Factors  not  taken  into  account  by  the  figures  of  merit, 

but  which  must  be  considered  in  selecting  a  material,  include 

acoustic  attenuation,  optical  transparency,  chemical  stability, 

mechanical  durability,  availability  of  an  established  technology 

for  producing  large  and  high  quality  boules ,  and  the  temperature 

coefficients  of  various  physical  constants,  especially  acoustic 
.  (23) 

velocity.  Optical  anisotropy  is  also  not  included  in  the 

(2k  25) 

figures  of  merit,  and  this  can  play  an  important  role.  ’ 

Table  1.1  below  provides  comparative  data  for  a  few 
selected  materials,  taken  from  reference  23.  The  measurement  of 
material  constants  and  the  identification  of  new  useful  materials 
is  an  important  area  of  acousto-optics  research. (26-29) 
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TABLE  1.1 

Selected  Acousto-Optic  Materials 


Material  and 
Propagation  Mode 

Useful 

Wavelengths 

(ym) 

Acoustic  Attenua¬ 
tion  at  500  MHz 
(dB/cm  •  GHz2) 

* 

M1 

# 

M2 

* 

M3 

Fused  Silica 

.2  -  2.5 

12 

1.0 

1.0 

1.0 

Water 

.2  -  .9 

2400 

4.6 

80.7 

18.52 

Dense  Flint  (SF-59) 

.46  -  2.5 

1200 

4.8 

12.2 

8.9 

Si:  Ltlll] 

1.5  -  10 

6.5 

25.6 

4.0- 

15.5 

Ge:  L[lll] 

2-20 

30 

1267. 

575. 

1370. 

GaP:  L [ 110 ] 

.6  -  10.0 

<  1 

74. 

29-5 

69. 

GaAs :  L[110] 

1-11 

30 

115. 

67. 

133 

LiNb03:  L[l00] 

.4  -  4.5 

.15 

8.3 

4.5 

7.5 

LiTaOy  L[001] 

.4  -  5 

.1 

1.4 

.8 

1.4 

a-HIO^:  L[001] 

.3  -  1.8 

10 

12.8 

55. 

31.1 

PbM  0,  :  L[001] 

0  4 

.42  -  5.5 

15 

15.5 

32.1 

26.0 

Te02:  L[001] 

.35-5 

15 

17.1 

22.1 

24.3 

Te02 :  S[110] 

.35  -  5 

290 

68 

793 

110 

* 

Normalized  to  fused  silica;  computed  from  bulk  wave  data. 
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1.1*  Acousto-Optic  Interaction  in  Bulk  and  Surface  Acoustic  Waves 

Until  recently  all  acousto-optic  devices  and  experiments 
have  employed  the  interaction  in  bulk  waves.  Because  of  the  con¬ 
siderable  research  and  development  effort  invested  in  bulk  wave 
device  technology  to  date,  it  is  to  be  expected  that  commercial 
applications  of  acousto-optics  will  be  dominated  by  bulk  wave 
devices  for  the  next  few  years  at  least.  This  section  discusses 
some  problem  areas  attendant  with  bulk  wave  devices  and  provides 
a  brief  introduction  and  review  of  surface  wave  acousto-optics. 

In  bulk  wave  devices  the  acoustic  wave  is  excited  by 
means  of  a  transducer  fabricated  from  a  plate  cut  from  a  piezo¬ 
electric  crystal,  across  which  the  driving  voltage  is  applied. 

A  variety  of  transducer  materials  are  available  for  this  purpose, 
the  most  widely  used  in  recent  practice  being  LiNbO^.^°^  Bonding 
of  the  transducer  to  the  acousto-optic  medium  is  the  most  impor¬ 
tant  technique  in  fabricating  a  bulk  wave  acousto-optic  device. 

The  use  of  an  organic  adhesive  such  as  epoxy  is  the 
simplest  bonding  technique.  Because  of  the  very  low  mechanical 
impedance  of  epoxy,  which  is  highly  mismatched  to  LiNbO^  and 
other  commonly  used  acousto-optic  materials,  the  thickness  of 
the  epoxy  layer  must  be  kept  to  a  small  fraction  of  an  acoustic 
wavelength.  Chang  has  obtained  layers  as  thin  as  0.1  ym 
which  have  given  good  results  up  to  frequencies  of  150  MHz. 

Above  this  frequency  other  techniques  are  required. 
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One  method  which  has  met  with  success  is  the  cold-weld 

( 31 ) 

bonding  technique  using  metallic  bonding  layers.  Indium 

is  the  most  commonly  used  metal  since  good  bonding  can  be  obtained 

using  low  bonding  pressure  ('v*  1000  psi)  for  a  few  seconds. 

However  indium  exhibits  a  rather  large  acoustic  attenuation  which 

precludes  its  use  at  high  frequencies.  Gold  provides  a  good 
(32) 

bond,  but  owing  to  its  extremely  high  mechanical  impedance 

it  is  mismatched  to  most  acousto-optical  materials.  Recently^0 ^ 

conversion  loss  of  less  than  1  dB  has  been  obtained  at  350  MHz 

using  aluminum  as  the  bonding  medium.  A  mechanical  impedance 

close  to  that  of  indium,  acoustic  attenuation  only  slightly 

higher  than  gold,  and  a  high  electrical  conductivity  combine 

favorably  in  this  material  for  use  at  high  frequencies. 

In  addition  to  good  transducer  bonding,  the  attainment 

of  high  acoustic  frequenices  with  bulk  wave  devices  requires  that 

the  transducer  plate  be  made  to  a  thickness  corresponding  to  the 

frequency  of  operation.  For  thicknesses  in  the  range  3-6  microns, 

( 32 ) 

mechanical  lapping  is  sufficient.  Sputter  etching  and  ion 
(33) 

milling  techniques  have  been  developed  to  further  reduce 

transducer  thickness,  .25  pm  having  been  obtained.  Stevenson 
(3M 

and  Hanak  have  reported  1  dB  conversion  loss  at  an  operating 
frequency  of  1  GHz  by  applying  ion  milling  and  indium-gold 
bonding  techniques. 

Another  problem  area  in  the  design  of  bulk  acousto¬ 
optic  devices,  and  surface  wave  devices  as  well,  is  that  of 
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designing  transducer  impedance  matching  networks.  In  the  hulk 

wave  case,  to  a  first  approximation,  the  transducer  acts  as  a 
( 35 ) 

large  capacitor.  Due  to  the  high  dielectric  constant  of 

LiNhO^j  the  impedance  level  of  transducers  for  most  devices  can 
be  on  the  order  of  a  few  ohms  for  frequencies  above  100  MHz. 

Thus  parasitic  inductances  can  dominate  the  transducer  impedance 
at  high  frequencies,  raising  the  loaded  Q  of  the  transducer  and 
narrowing  bandwidth.  The  use  of  multiple  series-connected  trans¬ 
ducers  is  often  necessary  to  reduce  transducer  capacitance  and 
raise  impedance  level.  Although  impedance  matching  networks 
based  on  filter  synthesis  have  been  developed,  broadband  matching 
is  often  obtained  by  resonating  the  transducer  with  a  shunt 
inductance  and  transforming  the  source  impedance  by  use  of  ferrite 
transformers  or  quarter  wave  transmission  lines. 

Recently  there  has  been  a  trend  toward  the  use  of  sur¬ 
face  acoustic  waves  for  acousto-optic  devices,  stimulated  by 

( 36 ) 

advances  in  surface  wave  microwave  acoustic  technology.  The 

motivation  for  the  development  of  surface  wave  acousto-optic 

devices  is  primarily  miniaturization,  ease  of  fabrication  using 

planar  technology,  and  compatability  with  the  emerging  integrated 
(37) 

optics  technology.  Further  advantages  of  the  use  of  surface 

( 38  i 

waves  over  bulk  waves  are  greater  transducer  flexibility  and 
the  reduction  in  one  dimension  of  acoustic  beam  diffraction  limi¬ 


tations. 


Theoretical  interest  in  surface  acoustic  waves  dates 
back  to  Lord  Rayleigh's  work  in  1885,  but  the  acousto-optic 
application  of  such  waves  was  not  practical  until  the  demonstra¬ 
tion  in  1988  of  an  efficient  means  of  generating  waves  on  the 
surface  of  the  acousto-optic  material  LiNbO^.^^  The  technique 
is  shown  in  Fig.  1.3.  Electro-mechanical  coupling  between  the 
interdigital  transducer  (IDT)  and  the  crystal  is  achieved  by 

means  of  the  piezoelectric  effect,  which  narrows  the  class  of 

(1*1) 

crystals  to  which  this  technique  can  be  applied.  The  use 

of  piezoelectric  materials  to  support  the  surface  wave  also  signi 

(1*2) 

ficantly  complicates  the  analysis  of  the  wave  propagation. 

In  the  figure  is  shown  a  single  transducer  with  two 
finger  pairs.  The  fingers  are  spaced,  from  center  to  center,  a 
distance  A/2  apart ,  ..where  A  is  the  acoustic  wavelength  at  the 
design  frequency.  For  YZ  LiNbO^,  operation  at  a  frequency  of 
100  MHz  corresponds  to  a  wavelength  of  approximately  35  pm. 
Generally,  the  electric-acoustic  conversion  efficiency  of  an 
interdigital  transducer  is  proportional  to  the  square  of  the 
number  of  finger  electrode  pairs.  However,  by  considering  the 
transit  time  for  an  acoustic  wave  to  travel  the  length  of  the 
transducer,  it  is  apparent  that  the  IDT  fractional  bandwidth  is 
given  by  the  reciprocal  of  the  number  of  finger  pairs.  In  YZ 
LiNbO^*  five-finger  pair  transducers  are  used  to  optimize  band¬ 
width  and  efficiency.  In  lieu  of  other  requirements,  the  trans¬ 
ducer  aperture,  length  of  finger  overlap,  is  determined  by 

(l*q) 

electrical  matching  considerations. 
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Fig.  1.3  Surface  Wave  Excitation 


fig  1.4  Interaction  Configurations 

(a)  Top  reflection  (b)  Top  transmission 

(e)  Back  reflection  (d)  Side  entry 
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A  variety  of  acousto-optic  applications  in  surface 

acoustic  waves  have  already  been  demonstrated.  These  have  included 

(1*1*) 

surface  wave  probing  for  microwave  acoustic  device  diagnosis, 

(1*5) 

real-time  convolution  using  direct  detection  and  optical 
heterodyning/^^  light  beam  modulation/**"^  image  scanning/*1*^ 
and  a  host  of  others.  Surface  acoustic  waves  have  also  found 
application  in  the  acousto-optic  investigation  of  layer  waves  /^^ 
Four  interaction  configurations  have  been  used  in  the 
acousto-optic  interaction  between  a  surface  acoustic  wave  and  an 
unguided  optical  wave,  shown  in  Fig.  l.U.  In  the  top  reflection 
configuration,  l.U(a),  light  is  reflected  from  the  surface  on 
which  the  acoustic  wave  propagates.  Diffraction  occurs  because 
of  the  slight  surface  aorrugation  associated  with  the  propagation 
of  the  acoustic  wave.  Although  the  interaction  is  weakest  in 
this  configuration,  it  is  a  useful  arrangement  for  the  non-contact 
probing  of  surface  waves. 

A  second  interaction  geometry  is  shown  in  l.U(b),  the 
top  transmission  configuration.  The  incident  laser  light  is 
directed  nearly  normal  to  the  acoustic  propagation  surface,  and 
the  diffracted  light  exits  through  the  bottom  of  the  crystal. 
Because  the  strain  components  of  the  surface  acoustic  wave  decay 
with  distance  from  the  crystal  surface,  this  configuration  has  a 
very  short  interaction  length,  on  the  order  of  an  acoustic  wave¬ 
length.  A  two- fold  improvement  in  the  interaction  length  is 
achieved  through  the  scheme  shown  in  l.l»(c),  the  back  reflection 
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configuration.  The  light  is  here  incident  as  in  the  top  transmission 

case,  but  after  exiting  from  the  back  side  the  light  is  reflected 

back  by  mirror  through  the  crystal  for  a  second  pass  through  the 

acoustic  wave.  The  periodic  corrugation  of  the  crystal  surface  due 

to  the  acoustic  wave  provides  a  slight  additional  phase  contribution 

to  that  due  to  the  internal  effect.  The  top  transmission  and  back 

(l+o  50) 

reflection  configurations  have  received  considerable  attention. 

The  side  entry  technique.  Fig.  l.U(d),  provides  a  vast 

improvement  in  interaction  length  over  the  other  configurations. 

The  incident  light  in  this  case  enters  from  the  side  of  the  crystal, 

near  the  crystal  top  surface,  and  propagates  parallel  to  the  surface 

on  which  the  acoustic  wave  is  propagating.  The  interaction  length 

for  the  side-entry  mode  of  operation  is  given,  therefore,  not  by  the 

depth  of  penetration  of  the  acoustic  wave  below  the  crystal  surface 

but  rather  by  the  acoustic  beam  width  which  is  determined  by  the 

aperture  of  the  interdigital  transducer.  For  a  transducer  on  YZ 

LiJTbO^  designed  to  match  to  50  fl  at  the  center  frequency,  this  is  an 

improvement  by  two "orders  of  magnitude  over  the  top  transmission  case. 

Recent  work  suggests  that  a  total  internal  reflection  condition  is 

achievable  with  the  side-entry  technique  which  provides  an  additional 

(1+7) 

enhancement  of  the  interaction  in  this  device  configuration. 

Although  optimum  optical  alignment  is  more  difficult  to  achieve  with 
the  side-entry  technique  than  with  the  other  configurations ,  it  is 
generally  quite  easy  to  exceed  the  interaction  strength  possible  with 
the  other  interactions  with  even  a  non-optimal  alignment. 


MHUMiH 
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Kramer  has  used  the  side-ent  ry  technique  for  mapping  the 
decay  of  the  Rayleight  wave  stran  fields  in  LiHbO^  and  quartz, 
and  has  reported  on  a  real-time  signal  processor  using  this  configura¬ 
tion.  More  recent  work  at  RPI  has  exploited  this  device  configura¬ 

tion  for  a  variety  of  applications . (52-56) 

One  feature  common  to  all  acousto-optic  devices  employing 
surface  acoustic  waves  is  that  the  figures  of  merit  are  larger  than 
those  for  hulk  devices  employing  the  same  material.  This  is  due  to 
the  fact  that  the  surface  wave  acoustic  velocity  is  invariably  smaller 
than  the  corresponding  bulk  wave  velocity,  on  the  order  of  a  factor  of 
1.5*  This  corresponds  to  a  factor  of  three  improvement  in  the  figure 
of  merit  M^.  Another  advantage  of  surface  acoustic  waves  over  bulk 
waves  is  the  very  small  acoustic  beam  waist  aspect  ratio,  allowing 
for  very  efficient  diffraction.  However,  of  the  four  interaction 
configurations  shown  in  Fig.  l.U,  only  the  side-ent  ry  technique  takes 
advantage  of  this  fact. 

The  investigation  of  the  acousto-optic  interaction 

between  surface  acoustic  waves  and  guided  optical  waves  is  an 

(51) 

especially  active  area  of  applied  optics.  Beam  deflection, 
mode  conversion, and  convolution^^  have  been  reported  to 
date.  Optical  beam  quality  and  in/out  coupling  efficiency  is 
good  for  small  optical  apertures,  which  makes  these  devices  well 
suited  to  beam  steering  applications,  especially  considering  the 
potential  low  power  requirements  of  these  devices.  Coupling  over 
a  large  optical  aperture  remains  a  serious  problem,  and  for  the 
present  this  places  a  severe  limitation  on  the  practical  application 
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of  these  devices  to  signal  processors  requiring  a  large  time- 
bandwidth  product,  since  this  is  proportional  to  the  optical 
aperture. 


1.5  Competitive  Technologies 

In  some  application  areas  alternative  methods  to  the 
acousto-optic  interaction  are  available.  In  the  case  of  con¬ 
trollable  light  deflection,  mechanical  and  electro-optic  means 
have  proven  useful.  The  mechanical  deflector,  a  mirror  or  prism, 
is  driven  by  piezoelectric,  galvanometric  or  magnetostrictive 
means.  Its  advantages  are  large  number  of  resol veable  spots, 
large  deflection  angle,  and  low  optical  loss.  Deflection 
accuracy  and  scanning  speed  are  problem  areas  which  have  yet 
to  be  overcome. 

Electro-optic  deflectors  inherently  exhibit  fast 
responses  and  have  been  extensively  investigated.  With  the 
exception  of  nitrobenzene  Kerr  cells  for  digital  light 
deflection,  the  unavailability  of  a  high  quality,  large  area 
electro-optical  material  having  the  capability  of  low  power 
operation  has  prevented  the  use  of  the  electro-optic  deflector. 

In  the  area  of  light  modulation  for  the  purpose  of 
optical  communications ,  the  direct  electronic  modulation  of  the 
LED  or  solid  state  laser  is  sufficient  and  practical  up  to  fre¬ 
quencies  of  about  150  MHz.  To  achieve  the  wide  bandwidth  poten¬ 
tial  of  single  mode  optical  fibers,  the  acousto-optic  interaction 
may  prove  useful  at  higher  frequencies. 
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For  purposes  of  signal  processing,  conventional  elec¬ 
tronic  and  real-time  digital  computing  methods  are  available  for 
low  frequency  applications.  In  certain  situations,  including 
processing  of  optical  images,  charged  coupled  devices  might  prove 
competitive.  At  higher  frequencies  where  the  acousto-optic 
interaction  might  usefully  be  applied,  other  techniques  have  been 
demonstrated  as  well.  Promising  among  these  are  those  devices 
employing  nonlinear  coupling  mechanisms  between  a  surface 

acoustic  wave  propagating  on  a  piezoelectric  substrate  and  a 

( ) 

semiconductor  placed  in  close  proximity.  Ultimately,  both 

device  performance  as  well  as  economic  factors  determine  the 
appropriateness  of  a  given  technology  for  a  given  application. 


CHAPTER  2 


THEORY  OF  THE  ACOUSTO-OPTIC  INTERACTION  IN  SURFACE  ACOUSTIC  WAVES 

2.1  Introduction 

The  theory  of  the  acousto-optic  interaction  in  surface 
acoustic  waves  is  developed  in  this  chapter  and  applied  to  the 
material  LiNbO^.  As  much  work  has  been  done  in  this  area  by  other 
authors,  an  overview  indicating  the  contributions  of  this  chapter 
will  be  helpful. 

Section  2.2  provides  the  foundation  for  the  subsequent 
analysis.  In  that  section  the  formulation  of  Pockels^^  is  presented. 
This  formulation  couples  a  mechanical  deformation  to  a  material’s 
optical  properties  by  relating  the  optical  permittivity  tensor  to 
elastic  strain.  Account  is  also  taken  of  the  fact  that  in  piezo¬ 
electric  materials,  strains  are  associated  with  electric  fields  which 
can  also  perturb  the  permittivity  tensor.  This  electro-optic  contri- 
bution  to  the  acousto-optic  effect  has  been  discussed  by  Coquin. 

In  Section  2.h  the  wave  equation  of  classical  electromagnetic 
theory  is  used  to  derive  a  set  of  coupled  partial  differential  equa¬ 
tions  to  describe  the  optical  field.  The  strain  and  electric  fields 
associated  with  a  material  deformation  appear  in  these  equations, 
introduced  through  the  permittivity  as  discussed.  Similar  coupled 
equations  have  been  given  by  Klein  and  Cook^"^  who  base  their  analysis 
on  a  scalar  wave  equation,  Chang  in  analyzing  the  case  of  a  bulk 
wave  in  a  non-piezoelectric  material,  and  Lean^"^  in  analyzing  the 
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surface  wave  case  by  means  of  a  Green's  function  technique.  In 
particular.  Lean  investigates  the  top  transmission  case  including 
a  contribution  from  the  surface  corrugation.  Using  the  acoustic 

ifQ  \ 

parameters  provided  by  Spaight  and  Koerber,  the  equations  are 

provided  for  the  specific  case  of  LiNbO^. 

In  Section  2.5  solutions  are  found  for  the  formulated 

problem  for  the  side-entry  configuration.  The  problem  had  not 

previously  been  analyzed  using  actual  Rayleigh  wave  parameters, 
(71) 

although  Kramer  provides  similar  results  by  considering  a 
surface  wave  in  an  acoustically  isotropic  material.  Numerical 
results  are  provided  in  Section  2.6  and  are  compared  with 
Kramer's  experimental  findings. 


The  change  in  a  material's  refractive  index  due  to  a 
deformation  is  generally  referred  to  as  the  photoelastic  effect. 

In  this  chapter  we  employ  a  modification  of  the  phenomenological 
formulation  introduced  by  Pockels  to  describe  this  effect.  In 
his  formulation  a  mechanical  strain  gives  rise  to  a  proportional 
perturbation  of  the  material's  relative  phase  tensor  B. 

The  phase  tensor  elements  are  defined  through  the 
linear  relation 

Eh  =  ~  B.j  Dj  ,  (i,j  =  1,2,3)  (2 

which  relates  the  components  of  the  electric  and  displacement 
fields  of  an  optical  wave  propagating  in  a  crystal.  Equation  (2.1) 
is  written  in  rationalized  MKS  units,  which  are  used  throughout, 
in  which  eQ  is  the  permittivity  of  free  space.  Throughout  this 
chapter,  unless  otherwise  noted,  the  Einstein  notation  is 
employed  in  which  summation  over  repeated  subscripted  indices  is 
implied.  According  to  the  formulation  of  Pockels,  the  change  in 
the  relative  phase  tensor  elements  at  optical  frequencies 


is  given  by 
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where  P. „  are  the  elements  of  the  fourth  rank  elasto-optic  tensor 
ljkk 

and  is  the  infinitesimal  strain  defined  below. 

If  u  represents  the  vector  displacement  of  a  material 
point  due  to  a  deformation,  and  are  the  magnitudes  of  the  com¬ 
ponents  in  an  orthogonal  coordinate  system  with  axes  x^  (i  =  1,2,3), 
then  in  the  so-called  "small  deformation  theory"  of  continuum  mechanics 
in  which  the  displacement  gradients  3u^/3x^  are  small  compared  to 
unity,  the  infinitesimal  strain  tensor  is  given  by 


1  3uJl  3uk 

Sk£  =  2  +  3^  *  (2,3) 

While  we  take  Eq.  (2.3)  as  the  definition  of  strain  for  our  pur¬ 
poses,  is  in  fact  the  infinitesimal  limit  to  which  both  the 
Lagrangian  and  Eulerian  strain  tensors  of  finite  deformation  theory 
reduce.  While  in  finite  deformation  theory  a  distinction  must  be 
made  as  to  whether  the  coordinate  system  is  fixed  in  space  or  fixed 
in  the  body  of  the  material  continuum,  no  such  distinction  need  be 
made  in  the  present  "small  deformation"  treatment. 

In  piezoelectric  materials,  in  which  the  strain  is  coupled 
via  the  piezoelectric  effect  to  an  associated  electric  field,  a 

deformation  can  also  give  rise  to  a  change  in  refractive  index  by 

f  63 ) 

means  of  the  electro-optic  effect.  It  was  pointed  out  by  Coquin v 
that  in  lithium  niobate,  which  is  highly  piezoelectric,  this  can 
be  a  large  effect,  as  important  as  the  direct  photoelastic  effect. 

Thus  we  must  write 


AB  =  P.jjj  S  p  +  rfa)  E  (a) 
ij  ijk£  k*.  ljm  m 


(2.M 


2b 


(  3  ) 

where  r. .  is  the  linear  electro-optic  tensor  element  measured  at 
ljm 

( E ) 

constant  strain,  is  the  elasto-optic  tensor  element  at  con¬ 

stant  electric  field,  and  E  is  the  mt'*1  component  of  the  elec¬ 
tric  field  associated  with  a  deformation  in  the  piezoelectric 

material.  The  superscript  (a),  referring  to  acoustic,  distinguishes 
(a) 

from  the  electric  field  associated  with  the  light,  which  is 
oscillatory  at  a  frequency  much  higher  than  that  associated  with 
piezoelectric  phenomena. 

Recently  Nelson  and  Lax^^  have  pointed  out  that  the 
phenomenological  formulation  of  Pockels  is  incomplete  and  neglects 
rotation  effects  arising  from  shear  waves  which  should  be  included 
in  strongly  birefringent  materials.  Thus  in  place  of  the  strain, 
a  symmetric  tensor,  the  displacement  gradient  which  is  asymmetric 
is  taken  by  them  as  the  more  basic  quantity.  Equivalently,  we 
could  include  a  contribution  to  the  optical  relative  inverse 
dielectric  tensor  (impermeability  tensor)  proportional  to  the 
rotation  tensor,  so  that 


AB. , 
ij 


=  P.  ,,  o  S  „  +  q  n  +  r,  .  E 
ljkJt  kZ  ijst  st  ljm  m 


(a) 


(2.5) 


where  the  rotation  tensor  U  is  the  antisymmetric  second  rank 
tensor  with  elements  defined  by 
,  3u  3u, 

o  s  i  ( — —  _  _ —) 

st  2  '  3x  3x  ' 

t  s 


Except  for  the  material  rutile. 


(66) 


(2.6) 


experimental  observation  has 


not  been  made  of  the  existence  of  nonvanishing  tensor  elements 
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^ijkJT  We  are  *nterested  a  numerical  solution  for  the  acousto¬ 
optic  effect  in  lithium  niobate,  we  do  not* include  the  rotation  con¬ 
tribution  in  subsequent  analysis,  although  such  inclusion  is  easily 
facilitated  within  the  framework  to  be  developed.  It  should  be 
specifically  noted,  too,  that  the  present  formulation  is  based  upon 
linear  piezoelectric  constitutive  relations,  the  linear  electro¬ 
magnetic  constitutive  relations,  and  linear  elasto-optic  and  electro¬ 
optic  relations. 

In  the  description  of  the  acousto-optic  interaction  to 
be  developed,  our  interest  is  not  in  how  the  impermeability  tensor 
is  perturbed  by  the  mechanical  deformation,  but  rather  how  the 
optical  dielectric  tensor  itself  is  perturbed.  Writing  the  rela¬ 
tive  permittivity  as 


=  e. ,  +  Ae.. 

ij  ij 


(2.7) 


where  is  the  unperturbed  tensor  element  and  Ae„  is  the  per¬ 
turbation,  it  follows  from 


D.  =  e.  .  E. 
i  ij  j 

and  Eqs.  (2.1)  and  (2.1+)  that,  to  linear  terms  in  S,  „  and  E 

kJi  m 

Ae..  =  -  (e.°.  )2  [P.,,„  S,.  +  r.,  E  . 

ij  li  ijk£  kl  ljm  m 


(2.8) 


(2.9) 


The  values  of  the  unperturbed  dielectric  tensor  elements  e^j  are 
those  appropriate  at  optical  frequencies.  In  the  coordinate  system 
in  which  the  unperturbed  dielectric  tensor  is  diagonal,  these 
elements  are  given  by 
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o  o  ,  r  (i )  i 

e.  .  =  e. ,  6 .  .  =  [n  ]  , 

iJ  iJ  ij 


(2.10) 


where  6..  is  the  Kronecker  delta,  and  is  the  refractive  index 

U 

for  light  polarized  along  the  principal  axis  x^.  In  lithium  nio- 
bate,  an  optically  uniaxial  crystal,  t  e33* 


2. 3  The  Acousto-Optic  Interaction 


As  discussed  in  Chapter  1,  when  an  ultrasonic  wave 
propagates  through  a  transparent  material,  the  refractive  index 
of  that  material  is  perturbed  by  the  presence  of  that  wave  via 
the  photoelastic  effect.  The  refractive  index  perturbation  is 
periodic  in  the  direction  of  acoustic  propagation.  Light  passing 
through  the  medium  can  undergo  diffraction  as  a  direct  result  of 
this  perturbation.  This  effect  is  referred  to  as  the  acousto¬ 
optic  interaction.  If  the  interaction  rotates  the  polarization 
of  the  incident  light,  we  refer  to  the  interaction  as  a  polari¬ 
zation  rotating,  or  simply  rotating,  interaction.  If  the  inter¬ 
action  does  not  rotate  the  light  polarizaiton,  we  refer  to  the 
interaction  as  nonrotating.  In  the  former  case  the  refractive 
index  for  the  incident  and  diffracted  light  will  generally  not  be 
equal  in  an  optically  anisotropic  material.  Such  would  be  the 
case  in  lithium  niobate  for  a  diffraction  process  coupling  light 
polarized  in  the  x^  direction  with  light  polarized  along  either 


of  the  other  axes.  In  both  the  rotating  and  nonrotating  cases, 
the  diffracted  light  is  generally  shifted  in  its  frequency  from 
the  incident  light  by  an  integral  multiple  of  the  acoustic 
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frequency.  For  this  reason  some  authors^  ’  view  the  inter¬ 


action  as  a  parametric  interaction  in  which  the  optical  wave  mixes 


with  the  acoustic  wave  to  generate  polarizion  waves  at  the  com¬ 


bination  frequencies  which  then  radiate  at  these  new  frequencies. 


Even  though  the  linear  constitutive  relations  and  linear  elasto- 


and  electro-optic  relations  are  employed,  the  acousto-optic 


interaction  is  inherently  nonlinear.  The  mixing  of  frequency 


occurs  in  our  present  theory  due  to  a  term  in  the  wave  equation 


for  the  optical  field  which  involves  the  product  of  the  optical 


field  with  the  perturbed  part  of  the  permittivity. 


2.1*  Coupled  Mode  Equations  and  Phase  Matching  Considerations 


Starting  with  Maxwell's  equations  the  following  wave 


equation  is  obtainable  which  describes  the  electric  field  of  the 


optical  wave  in  the  crystal.  The  details  of  this  derivation  are 


shown  in  the  Appendix. 


-  V2  E.  +  [V(V  •  E)].  =  -  V  -2-=-  e..  E, 

i  i  o  9t2  ij  j 


(2.11) 


It  is  further  shown  in  the  Appendix  that  the  gradient  term  in  Eq. 


(2.11)  is  small  compared  to  all  other  terms  and  can  be  neglected. 


The  resultant  equation,  which  provides  the  basis  for  the  present 


investigation  of  the  acousto-optic  interaction,  relates  the  optical 


field  to  the  material's  permittivity  which  is  perturbed  by  the 


acoustic  wave.  Written  in  a  cartesian  coordinate  system,  employing 


I 


indicial  notation 


where  is  the  i  component  of  the  electric  field  associated 


with  the  light  propagating  in  the  material,  and  e..  is  an  element 

^  J 


of  the  permittivity  tensor,  given  explicitly  by  Eq.  (2.?). 

With  reference  to  the  geometry  of  Fig.  2.1,  the  acoustic 
strain  and  associated  electric  field  can  be  written  as 


Sk£  =  2  +  skZ  ^ 


and 


<*>  -  i<.  ♦  * 


m 


(e  +  e  ) 
2  m  m 


in  which  the  asterisk  *  denotes  complex  conjugate,  and 

Ja^y)  jUat-kaz) 


skA  =  Fk£  (y)  e 


e  =  G(y)  e  “ 
m  m 


jtra(y)  j(u)  t-k  z) 


a  a 


( 


where  j ,  Gffi,  a^.  and  b^  are  real-valued  functions  of  the  vari¬ 
able  y=  x^.  The  subscript  "a"  on  quantities  in  Eqs.  (2.15)  and 
(2.1 6)  refers  to  the  acoustic  wave.  The  coordinates  (x,y,z)  and 
(x^jX^jX^)  are  used  interchangeably  throughout  the  chapter.  While 
the  latter  enable  us  to  use  the  compact  indicial  notation,  the 
former  are  conventionally  used  in  specifying  crystallographic  cuts 
The  imaginary  unit,  j  =  i£T,  is  easily  distinguished  from  a  sub- 
or  superscript  j,  used  as  an  index,  from  its  context. 
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While  we  are  concerned  with  the  solution  to  Eq.  (2.12), 
subject  to  boundary  conditions,  for  practical  cases  in  which  the 
width  W  of  the  incident  light  beam  is  of  finite  extent  in  the  z- 
direction,  or  the  acoustic  disturbance  itself  is  of  finite  spatial 
extent  on  the  crystal  or  is  not  of  constant  amplitude,  such  effects 
will  not  be  taken  into  explicit  account  in  the  following  analysis. 
The  effect  of  a  finite  optical  aperture  is  well  understood  and  is 
discussed  in  optics  texts  treating  diffraction.^^  In  subse¬ 
quent  chapters,  in  which  applications  of  the  interaction  are  dis¬ 
cussed,  the  effect  of  a  finite  optical  aperture  is  introduced  where 
required.  The  case  in  which  the  acoustic  disturbance  is  not  of 
constant  amplitude  or  is  of  finite  spatial  extent  on  the  crystal 
is  of  crucial  importance  to  the  operation  of  all  the  devices 
described.  The  modifications  required  to  apply  the  results  of 
this  chapter  in  describing  those  devices  are  explicitly  noted 


for  each  device  application  in  the  appropriate  section. 

It  follows  from  Eqs.  (2.9),  (2.13),  (2.lU),  (2.15), 
and  (2.16)  that  the  perturbed  relative  permittivity  is  of  the  form 


-  Vy) 


J  ( 0)  t  -  k  z) 

a,  a 


(2. IT) 


in  which  and  h^  are  real  valued  functions  of  y.  It  is  shown 
in  the  next  section  that  for  the  acousto-optic  interaction  with 


surface  waves  in  lithium  niobate,  only  the  functions  ^"33 


are  of  practical  interest  for  the  geometry  shown  in  Fig.  2.1. 
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For  a  z-propagating  Rayleigh  wave  on  a  stress  free  y-cut 


lithium  niobate  substrate,  Ae^  has  the  analytic  form 


l  o  2  (1)  ~iaiKy 

V  =  (V  Baka  tCij  e  +  V  e 


(2)  _“Ja2kay 


+  C 


(3)  ~JOt3kay 
ij  6 


]  +  C.C. 


(2.18) 


where  B&  is  the  amplitude  of  the  acoustic  wave,  and  the  other  con¬ 
stants  characterize  the  material.  The  coefficients  are  the 

ij 


(m)s 

ij 


sum  of  an  elasto-optic  and  an  electro-optic  contribution,  C 
and  respectively.  All  constants  are  provided  in  Section 

2.5  of  this  chapter. 

As  a  solution  to  the  wave  equation  (2.12)  we  take  a  plane 
wave  expansion  of  the  form 


E.(x,y,z,t)  =  Z  (x,y)  e  n  n 

Y\  =  —00  n 


j(w  t-k^  •  r) 


(2.19) 


n  = 


in  which  is  the  amplitude  of  the  n^h  diffracted  order,  "r  is 

the  position  vector,  u>n  the  frequency  of  the  ntk  order,  and  ^ 
the  wave  vector  associated  with  the  n  1  order  optical  mode.  The 
superscripts  on  the  mode  amplitude  and  wave  vector  indicate  that 
these  depend  on  the  particular  electric  field  component  (i.e. 
light  polarization).  The  frequency  u>a  and  wave  vector  '  obey 


b> 


=  0)  ,  +  U) 

n-1  a 


(2.20) 


*  (68) 

Based  on  Spaight  and  Koerber. 


in  which  u)  is  the  frequency  of  the  incident  light  and  k  is  the 

O  8i 

acoustic  wave  vector.  The  wavevector  of  the  incident  light  is 
«•*  /  \ 

denoted  by  kQ  p  ,  where  p  is  either  2  or  3.  It  follows  from  Eq. 
(2.20)  that 

ui  =  to  +  n  U)  (; 

n  o  a 


Referring  to  the  geometry  of  Fig.  2.1, 


(i) 


cos 


(i) 


n 


x  +  k 


(i) 


sin 


u) 


(; 


where  all  wavevectors  and  angles  in  Eq.  (2.23)  are  measured  in  the 
medium.  Using  Eqs.  (2.2l)  and  (2.23),  conservation  of  momentum 
in  the  x-direction  provides 

k„u>  e„(i)  -  •‘Ji’  •  « 

and  conservation  in  the  z-direction  yields 


k  o-in  fl  ^ 

k  sin  u 

n  n 


=  k 


(j) 


n-1 


sin  6 


(j) 


n-1 


+  k 


C 


In  terms  of  the  incident  light  polarized  along  x  axis 

P 


(p  =  2,3),  Eqs.  (2.24)  and  (2.25)  can  be  written 
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k  ^  sin  0  ^  =  k  ^  sin  9  ^  +  n  k 

n  n  o  o  a 


(2.27) 


The  angles  can  be  related  to  quantities  in  free  space  by  means  of 
Snell's  Law,. 


.  _  (i )  .  A  (i) 

sin  0  _  -  n  sin  o 

nf  n 


(2.28) 


in  which  the  subscript  f  denotes  free  space.  Similarly,  wave 
vector  magnitudes  can  be  related  to  free  space  quantities  by 


n 


(i) 


knf  = 


(2.29) 


in  which  the  wave  vector  magnitude  for  free  space  propagation  k^, 
is  related  to  the  free  space  wavelength  A  by 


k 


nf 


(2.30) 


( i ) 

It  should  be  observed  that,  in  general,  n  depends  upon  the 

direction  of  the  propagating  wave.  We  assume  that  all  angles 

considered  are  sufficiently  close  that  the  refractive  index  may 

* 

be  treated  as  a  constant  for  each  polarization.  For  incident 
light  not  polarized  along  the  y-  or  z-axis,  the  light  must  first 
be  decomposed  into  components  along  those  axes.  Each  component 
is  then  treated  separately,  the  results  being  combined  at  the 
exit  face  of  the  crystal. 

With  the  definition 


,(i)  (x,y)  =  *  ' (i)(x,y)  e 


-3  K 


(i) 


cos 


(i) 


(2.31) 


The  assumption  is  good  for  the  experiments  reported  here,  but  becomes 
invalid  at  higher  acoustic  frequencies. 


the  field  can  be  written 


E. 

1 


00  /.»  j  [w  t-  (k  ^  sin  0  ^  +  n  k  )z] 

-  .  (i),  v  °  n  o  o  a 

Z  ij>n'  (x,y)  e 


(2 


Substituting  Eq.  (2.32)  into  the  wave  equation  (2.12)  and  equating 
terms  at  the  light  frequency  (oi  +  n  to  )  provides  a  coupled  system 

O  0. 

of  linear  second  order  coupled-mode  partial  differential  equations. 


-4  *  (i>  ♦  -4  *  (i)  ♦  Ar  e.°.  -  Oc  (1)  sin  8  (i)  *  „  k  n  * 

-  2  Tn  2  Tn  2  n  o  o  a  Yn 

3x  3y  c 


1  ^n  R  r  *^hii  .  (i)  ~^hii  ,  (i), 

=  Bii[*  *n-l  *e  Vl  1 

c 


i  (0 

1  n 


[eJh«  *  <j>  ♦  9‘3hiJ  *  (J> 

n-1 


2  c2  ij 


n+1 


(2 


where  i  ^  j  in  this  equation.  Use  has  been  made  of  the  fact  that 
in  the  coordinate  system  employed,  the  unperturbed  permittivity 
tensor  is  diagonal. 

Several  comments  are  in  order  regarding  Eq.  (2.33). 

It  must  be  pointed  out  first  that  this  represents  two  equations 
which  are  coupled,  since  i  =  2,3.  The  term  corresponding  to  i  =  1 
is  neglected  since  the  case  with  which  we  shall  be  concerned,  all 
light  propagation  is  along  the  x-axis,  or  nearly  so,  with  the 
result  E^  0*  Further,  for  the  particular  surface  wave  considered 
a  z-propagating  Rayleigh  wave  on  a  y-cut  crystal,  there  is  no  x- 
component  of  mechanical  strain,  thus  j  =2,3  only.  The  diffracted 


orders  (modes)  are  coupled  to  adjacent  orders  through  the 

subscript.  This  coupling  can  be  to  light  of  both  the  same  and 

also  orthogonal  polarization,  as  provided  by  the  superscript. 

The  first  terra  on  the  right  hand  side  provides  for  the  case  in 

which  the  adjacent  modes  are  of  the  same  polarization.  The  second 

term  on  the  right  hand  side  provides  for  coupling  between  different 

polarizations  and  allows  for  the  possibility  that  the  refractive 

indices  for  the  adjacent  modes  can  be  different.  In  Eq.  (2.33), 

coupling  is  only  between  adjacent  orders,  i.e.  <J>  is  coupled 

to  and  This  arises  from  the  fact  that  in  Eq.  (2. IT) 

n+1  n-1 

has  a  simple  sinusoidal  time  dependence.  In  the  more 
general  case  that  Ae.^  is  periodic  in  time,  expressed  as  a  Fourier 
sum  of  sinusoids,  coupling  would  be  possible  to  other  than  just 
the  adjacent  mode. 

Before  proceeding  to  construct  solutions  to  Eq.  (2.33) 
subject  to  boundary  conditions,  considerable  physical  insight  as 
well  as  an  important  mathematical  simplification  can  be  achieved 
by  a  more  careful  consideration  of  the  phase  expressions,  Eqs. 

(2.23)  through  (2.27).  We  prove  first  that  in  all  practical 
cases  there  cannot  exist  polarization  rotation  between  axes  of 
different  refractive  index  in  any  multiple  order  diffraction, 
such  as  Raman-Nath  diffraction.  We  prove  this  by  assuming  such 
an  interaction  can  occur  and  show  that  this  implies  an  incon¬ 
sistency. 

Consider  coupling  between  the  incident  light  with 
polarization  along  the  x^  axis  and  the  first  order  with  polarization 


along  x^.  Squaring  Eq.  (2.26)  and  substituting  from  Eq.  (2.27) 
provides  for  n  =  1, 


sin  9 


(i)  _  ka  +  °f  ,  (i)2  (p) 


2n(1>  k 


(2.31*) 


In  terms  of  the  acoustic  frequency,  Eq.  (2.3l*)  can  be  written 


sin  9. 


X  2 

o,  v  , .  k2  ,  s2 

uy— tf. +  — f— <■>  -"(p) 

n  v  X  f 

a  Of  a 


(2.35) 


in  which  f  and  v  are,  respectively,  the  acoustic  frequency  and 

Si  EL 

velocity.  Equation  (2.35)  is  a  valid  and  useful  result  for  an 
acousto-optic  interaction  involving  only  two  modes,  such  as  Bragg 
diffraction  which  involves  an  incident  and  a  single  diffraction 
mode.  It  is  a  feature  of  Eq.  (2.35)  that  there  is  a  critical 
frequency  f'  at  which  the  two  terms  in  brackets  on  the  right  hand 
side  are  equal.  This  frequency  is  given  by 


(2.36) 


For  f  <  f ' ,  the  second  term  in  the  brackets  of  Eq.  (2.35)  is  the 

EL 

larger  of  the  two  terms,  in  magnitude.  This  is  generally  the  case  for 

a  polarization-rotating  interaction  involving  a  surface  acoustic 

o 

wave  in  LiNbO_.  For  X  =  6328  A,  f'  =3-5  GHz.  Experiments  reported 
^  °f 

here  were  well  below  this  frequency. 
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Having  established  the  foregoing  it  is  easily  shown  that 
multiple  diffractions  involving  polarization  rotation  cannot  occur 
below  1/2  f ' .  If  it  could  occur  then  coupling  to  the  seccr.d  order 
would  provide 


sin  9, 


(p)  _ 


[f  +  - 

v  k  X  ^  f 

a  of  a 


(2.37) 


which  is  derived  in  exactly  the  same  manner  that  Eq.  (2.35)  was. 

Below  frequencies  1/2  f'  the  sign  of  sin  0  is  given  by  the 
2  2 

sign  of  [n^  -  n^  ],  while  the  sign  of  sin  0  is  of  the 

opposite  sign.  If  multiple  order  diffraction  involving  polariza¬ 
tion  rotation  could  occur,  both  angles  would  be  of  the  same  sign. 
For  this  reason,  in  studying  Raman-Nath  diffraction  in  highly 
birefringent  materials,  we  need  not  include  the  term  involving 
R^j  (i  i  j)  in  Eq.  (2.33). 

For  an  interaction  involving  only  a  single  diffracted 
order,  such  as  Bragg  diffraction,  polarization  rotation  is  not 
precluded  by  the  above  argument.  We  can  identify  two  possible 
cases.  In  the  event  that  the  diffracted  order  is  of  the  same 
polarization  as  the  incident  light,  Eq.  (2.33)  applies  with  the 
R^j  term  excluded.  In  the  event  that  the  diffracted  order  and 
the  incident  order  are  of  orthogonal  polar izaiton,  Eq.  (2.33) 
applies  with  the  R^  term  excluded.  The  case  that  the  diffracted 
order  is  a  superposition  of  both  polarizations  is  precluded  by 
Eq.  (2.35),  except  in  the  collinear  case  described  below. 
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It  follows  from  Eq.  (2.25)  that  for  the  case  of  an 
acousto-optic  diffraction  involving  only  one  diffraction  order, 
there  is  a  minimum  acoustic  frequency  at  which  the  interaction 
can  occur  if  the  diffracted  and  incident  light  are  of  orthogonal 
polarization.  This  frequency  occurs  for  the  case  in  which 
sin  =  sin  0  which  corresponds  to  the  incident  and 

diffracted  light  propagating  collinearly  with  the  acoustic  wave. 

This  minimum  acoustic  frequency  is  given  by 

f  •  =  y5-  (n(l)  -  n(p))  .  (2.38) 

min  A 

°f 

As  a  consequence  of  Eq.  (2.38),  which  assumes  that  the  incident 

and  diffracted  light  both  propagate  in  the  same  direction  as  the 

acoustic  wave,  if  n^  t  n^  then  the  interaction  can  only  occur 

if  n^^  >  n^.  That  is,  the  incident  polarization  must  be  along 

the  axis  of  smaller  refractive  index.  It  is  to  be  noted,  too, 

that  if  n^^  =  n^P\  such  as  for  an  optically  isotropic  crystal, 

or  for  light  propagation  along  the  c-axis  of  a  uniaxial  crystal, 

f  .  is  zero.  In  these  circumstances  the  interaction  is  possible 
mm  c 

for  acoustic  frequencies  down  to  the  dc  limit.  Finally,  had  we 
considered  the  case  in  which  the  acoustic  and  light  waves  propa¬ 
gate  collinearly  but  in  opposite  directions,  the  condition 
n^  >  n^  would  be  replaced  by  n^  >  n^. 

Equation  (2.19)  expresses  the  optical  field  in  the 
medium  as  a  linear  superposition  of  plane  waves  at  frequencies 

0)  =  u)  +  n  go  ,  (n  =  0,  +  1,  +  2,  ...).  For  light  polarized 
no  6L 


in  the  i"^  direction,  (i  =  2,3),  the  local  intensity  of  the 
order  is 


(i) 


=  k 


U)i 


(2.39) 


The  z-component  of  the  wavevector  for  this  wave  is  (k  sin  9  + 

o  o 

n  kj,  thus  the  sine  of  0  the  angle  between  the  x-axis  and 

a  n 

the  direction  of  propagation  of  the  wave  at  frequency  (on  is 


sin  6 

n 


(i) 


c(l^  [k  ^l)  sin  6  +  n  k  ] 

_ o _ o _ a_ 

<o  +  n  to 
o  a 


(2. bo) 


in  which  c^^  is  the  speed  of  light  in  the  medium  for  polarization 
along  the  x..  axis.  With  Snell's  Law,  and  using  the  approximation 
that 


n  to  «  co 
a  o 


it  immediately  follows  that 


n  X 


sin  6  =  sin  0  +  — t— 

nf  °f  A 


(2. Ill) 


(2.b2) 


in  which  the  subscript  f  refers  to  free  space.  The  quantities  0^, 

X  ,  and  A  are,  respectively,  the  angle  from  the  x-axis  of  the 
f 

direction  of  propagation  for  the  n  mode,  the  free  space  wave¬ 
length  of  the  incident  light,  and  the  acoustic  wavelength  which  is 
related  to  the  magnitude  of  the  acoustic  wavevector  by 


(2.U3) 


Of  special  interest  are  the  cases  corresponding  to 
Raman-Nath  and  Bragg  scattering.  In  the  former,  near  normal 


1*0 


incidence,  6  ^  0, 

°f  " 


n  X 


sin  6  'N*  — t — 

nf  —  A 

In  the  latter,  incidence  is  at  the  Bragg  angle,  X  =  0_,  where 

°f  B 

0D  is  defined  by 

.D 


(2.UU) 


•  1  °f 

sin  eB  -  '  2  ~r 


(2.1*5) 


assuming  the  refractive  indices  for  the  incident  and  diffracted 

wave  are  the  same.  Thus  the  diffracted  order  occurs  at  0  given 

1f 

by 


sin  0, 


2A 


(2.1*6) 


which  agrees  with  the  previous  result,  Eq.  ( 2 . 35 ) »  when  n 
,  as  expected. 


(i)  . 


2. 5  Simplification  of  the  Coupled  Mode  Equations:  Solutions  for 
Raman-Nath  and  Bragg  Diffraction 

In  the  preceding  section  the  coupled  mode  equations  were 
formulated,  and  arguments  based  on  the  conservation  of  pseudomomentum 
were  given  to  simplify  those  equations.  The  result  of  those  argu¬ 
ments  is  that  Eq.  (2.33)  becomes 


a2  .  (i)  .  a2 


.  u;  ,  a _  .  (i)  r  "n  o  ,,  .  «  li^  .  ,  \2-,  , 

2tn  +  -J  <*>n  +  t~  hi  ~  Sln  9o  +  n  ka)  ] 

9x  9y  c 


u 


(i)  a  (i) 


(i) 


1  “n  .  .  *Jhik(y)  <k)  t  -Jhik(y)  ,(k), 

*n-l  +  e  Vl1 


(2.1.7) 


1*1 


in  which  i  and  k  may  take  on  the  values  2  and  3.  If  the  diffracted 
light  is  of  the  same  polarization  as  the  incident  light,  which  is 
always  the  case  in  Raman-Nath  scattering,  i  =  k.  In  Bragg 
scattering,  if  the  incident  light  and  diffracted  light  are  of 
orthogonal  polarization,  i  /  k.  When  i  ^  k,  Eq.  (2.1*7)  repre¬ 
sents  a  coupled  set  of  partial  differential  equations.  When  i  =  k 
the  equations  decouple. 

As  was  mentioned  earlier,  the  field  E  associated  with 
the  incident  light  can  always  be  decomposed  into  components  along 
the  principal  axes  of  the  unperturbed  crystal.  Each  component 
can  then  undergo  acousto-optic  scattering  via  Eq.  (2.1*7),  and 
the  results  for  the  separate  incident  polarization  components 
linearly  superimpose.  It  is  therefore  only  necessary  to  consider 
the  case  in  which  the  incident  light  is  polarized  along  the  xp 
axis,  where  p  is  2  or  3.  Designating  by  q  the  orthogonal  axis, 
the  boundary  conditions  we  impose  at  x  =  0  are 

(p) 


(o,y)  =  D 


(2.1*8a) 


(p) 


m 


(o,y)  =  0  ,  (m  =  +1,  +2,  +3, . . . ) 


(2.  l*8b ) 


#n(q)  (o,y)  =  0  ,  (n  =  0,  +1,  +2,...)  (2.1*8c) 


where  a  uniform  plane  plane  with  amplitude  D  for  the  p-polarized 
component  has  been  assumed. 

We  proceed  now  to  further  simplify  Eq.  (2.1*7)  before 


separately  investigating  the  special  cases  of  Raman-Nath  and  Bragg 
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diffraction.  The  equation,  although  linear,  has  nonconstant 
coefficients  and  consequently  is  not  separable  in  the  variables 
x  and  y.  However,  because  no  mechanism  for  loss  has  been  intro¬ 
duced  in  our  analysis,  which  could  give  rise  to  spatial  amplitude 
modulation  of  the  light,  the  predominant  effect  of  the  y-derivative 
is  to  introduce  a  phase  variation  in  the  y-direction.  The  effect 
of  such  phase  variation  would  be  to  introduce  spreading  in  the 
y-direction  of  the  output  light  beam.  It  should  be  pointed  out 
that  in  all  applications  reported  in  the  following  chapter,  all 
the  light  in  a  given  diffraction  order  is  collected  by  the  photo¬ 
detector.  Any  variation  in  the  y-direction  is  not  important; 
only  the  total  integral  of  the  light  intensity  taken  over  the 
face  of  the  detector  is  of  practical  interest.  Because  of  this 
integration,  only  the  average  value  of  the  acoustic  wave  ampli¬ 
tude  taken  in  the  region  over  which  the  light  interacts  is  signi¬ 
ficant  in  determining  the  measured  light  diffracted  into  each 
order.  In  any  application,  what  is  measured  by  the  photodetector 
is  an  effective  value  of  the  permittivity's  perturbation  for  some 
interval  Ay. 


We  now  obtain  solutions  to  Eq.  (2.47)  subject  to  the 
boundary  constraints  (2.48).  We  investigate  first  the  case  in 
which  the  polarization  of  the  light  is  not  rotated  as  a  result 


1+3 


of  the  acousto-optic  interaction,  i  =  k.  Defining  Fn^  (x,y) 


by 


FnU)(x,y)  =  e 


"J  n  2  -Jnh..(y) 


e  11'"’  $nU)(x,y), 


(2.1*9) 


employing  Eq.  (2.31)  and  neglecting  the  y-derivative,  Eq.  (2.1+7) 
is  transformed  into  an  ordinary  differential  equation  in  which  y 
appears  as  a  parameter. 

2  j  k  cos  0  (l)  - - +  [  (k  cos  6  (i))2  -  e  ° 

o  o  dx  o  o  2  ii 

c 


+  (k  ^  sin  0  ^+nk  )2]  F  ^ 


<o 


-  i  i  ft. .  r  (<)  rp  (i)  p  ( i ) ■)  +  d  p  ( i ) 

-  J  2  _2  Rii(j;  n+1  "  Fn-1  J  +  ~2  Fn 

dx 


(2.50) 


With  the  change  of  variable 

co  2  1 

ft  =  £  Rii^  tk0(l)  cos  X 

c 


(2.51) 


Eq.  (2.50)  becomes 
(i) 


d  M 


dn 

(i) 


+  [M  ^  -  M  =  j  L  ^ 

n-1  n+1  u  n  n 

(i) 


in  which  M  and  L  are  defined  by 
n  n 


(2.52) 


Mn(l)(n)  =  Fn(i)  (x,y ) 


(2.53) 


and 


uu 


L  (1)  =  5-Vt  A  !(k  (I)  =os  6  (i))2  -  \  c° 
n  R. .  (y)  2  o  o  2  n 

11 w  w  c 

n 


2 

+  (k  ^  sin  8  ^  +  n  k  )  ] 


(2.5U) 


To  facilitate  comparison  of  our  results  with  those  of 
other  authors  we  define  and  8^  by 

/. x  k  2  L 

Qu)  -  _a _ 

9  „(1»k 


,u)  . 


n^  k  sin  9  ^ 1  ^ 

Of 


in  which  L  is  the  interaction  length  in  the  medium,  0  <_  x  _<  L, 
as  shown  in  Fig.  2.1.  Using  Eq.  (2.1*1)  we  can  write 


(2.55) 


(2.56) 


t,  (i>  *  aalii - 

n  -  R..(y)  k 


[(n  -  2  8V1;)] 


(2. 57) 


Klein  and  Cook,  in  considering  the  acousto-optic 
interaction  in  an  isotropic  medium  supporting  bulk  acoustic  waves, 
show  that  the  Raman-Nath  and  Bragg  regimes  correspond  to  the 
extreme  cases  «  1  and  >>  1,  respectively.  Similar 

schemes  for  differentiating  between  the  two  cases  have  been  proposed 
by  Willard^^)  Extermann  and  Wannier.^0^  Klein  and 


Cook  investigate  these  cases  analytically  and  provide 


numerical  solutions  for  the  transition  region,  Q  ^  1.  For  the 

the 

(12 


case  of  a  bulk  wave,  R.„(y)  =  constant,  Eq.  (2.52)  reduces  to  the 


set  of  coupled- mode  equations  first  derived  by  Raman  and  Nath. 


In  Raman-Nath  diffraction,  for  even  modest  non-normal 


,(i) 


angles  of  incidence,  can  be  very  large,  with  the  result  that 


(i) 


L  can  be  written  as 
n 


.  (i)  ,,  j.  n(1)  n  a(i)  e(1) 


In  terms  of  the  quantities  k  ,  k  ,  and  sin  0  ,  this  is  given  by 

a  Of  of 


U  n  sin  0  k 


,  U)  .  _ _ 

n  -  Ri.(y)  kQ 


Of  a 


f 


Substituting  Eq.  (2.59)  into  (2.52),  and  using  the  two  Bessel 
function  identities 


2  J  (x)  =  J  (x)  -  J  (x) 

n  n  -  1  n  + 1 


and 


Jn- 1  (x)  +  Jn  +  1  (x) 


—  J  (x) 
x  n 


it  follows  that 


Mn(l)  (n)  =  A  exp  [J  (  )f  n)  Jn  lsln  e  /  •!«  <k  -KT(y 

11*'  A  D  A  It* 


sin  0  k 

a 


Rii<y)  ko. 


sin  0  k 

i 


11  ’  o 


Of  a  of  li 


in  which  the  constant  A  is  determined  from  Eq.  (2.U8) 


r 
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In  the  limiting  case  of  a  bulk  wave,  R^y)  =  constant,  the  result 


is  of  the  same  for  as  given  by  Klein  and  Cook. 

The  intensity  of  the  n^*1  diffraction  order  is 

R. . (y)  k  sin  0  k 

!<*>.«.  D2J.2  [  \ . .  sin  (-— 

n  ip  n  sm  0^  k&  kQ  R..,.  (y) 


n)3 


(2-64) 


At  x  =  L,  this  reduces  in  the  case  of  normal  incidence  to 


InU)  =  6  in  ^  Jn2  C2  Rii(y) 
n  ip  n  d  ii 


LI 


(2.65) 


Using  the  expansion  for  integral  order  Bessel  functions  of  small 
argument , 


j  (a)  *  ir  <f)m  >  a  «  1  ,  (2.66) 

m  m!  2 

the  first  order  intensity,  normalized  to  the  incident  light,  is 
given  by 


[Ru(y) 


(2.67) 


A  result  analogous  to  Eq.  (2.67)  has  previously  been  given  by 
(71) 

Kramer,  arrived  at  by  analogy  to  the  Raman-Nath  diffraction 
in  bulk  waves  using  a  physical  argument.  His  theory,  however, 
assumed  an  acoustically  isotropic  material  and  did  not  take  the 
electro-optic  contribution  into  explicit  account.  . 

In  any  physical  experiment  the  photodetector  gives  a 


signal  proportional  to  the  integral  of  the  light  intensity  taken 


over  the  detector's  active  area.  The  geometry  for  an  idealized 


acousto-optic  experiment  is  shown  in  Fig.  2.2.  The  laser  beam 
of  width  w  in  the  y-direction  interacts  in  the  lithium  niobate 
a  depth  h  from  the  surface.  All  light  in  a  given  order  is 
collected  by  a  photodiode  after  focusing  by  a  lens.  Thus  the 
photodiode  output  is  given  by 


out 


=  C 


rect 


(i) 


(L,y)  dy 


(2 


in  which  a  uniform  plane  wave  of  width  w  has  been  assumed,  and 
is  a  constant  of  proportionality  which  characterizes  the 
detector,  given  by  the  product  of  the  photoelectric  conversion 
factor,  subsequent  amplification,  and  load  resistance.  The 
sampling  function  is  defined  by 


rect(x) 


(2 


Effects  such  as  boundary  reflections ,  cone-angle  due  to  lens 

focusing,  and  Gaussian  laser  beam  profile  will  modify  Eq.  (2.68). 

Cone-angle  and  Gaussian  profile  characteristics  can  be  taken  into 

account  by  properly  defining  h  and  w  and  have  been  treated  in 
(71) 

previous  work.  Boundary  reflections  can  be  eliminated  by 

proper  index  matching.  From  Eqs.  (2.68)  and  (2.69),  light  in 
the  first  diffracted  order  will  produce  the  following  photo¬ 


diode  signal: 


1*8  ! 
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_1_  r  f 
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rect  ( 


h  - 


WZ)  [Rpp(y)]2  dy 


(2.70) 


Following  the  method  of  Klein  and  Cook,  we  now  investi¬ 
gate  Eq.  (2.5 2)  for  the  case  of  Bragg  diffraction  in  which  the 
incident  and  diffracted  beam  are  of  the  same  polarization.  Unlike 
the  Raman-Nath  case  we  now  seek  a  solution  for  which  only  one 
diffraction  order  is  nonvanishing.  For  light  to  appear  in  some 
single  order  other  than  the  zero,  that  order  must  be  synchronous 
with  the  zero  order.  Thus  the  right  hand  side  of  Eq.  (2.52)  must 
have  the  same  value  for  the  diffracted  order  as  for  the  zero  order. 
Equation  (2.57)  provides  that  the  right  hand  side  is  zero  for  this 
order.  If  the  first  order,  n  =  1,  is  to  be  synchronous  with  the 
zero  order,  we  must  insist  that 


,(i)  =  1 


(2.71) 


For  n  =  -1,  we  must  have  8^  =  -1/2;  for  a  synchronous  second 


,  e(i)  =  1. 


A  condition  can  be  imposed  upon  &  * ^  to  make 


order 

L  ^  vanish  for  any  desired  order.  It  should  be  recalled,  how- 
n 

ever,  that  only  adjacent  modes  are  coupled  by  the  differential 
equation.  Thus  8^  =  +1  is  unacceptable.  $  ^  and  ^  are 
not  directly  coupled.  Moreover,  the  order  through  which  the  zero 
and  second  orders  are  coupled,  the  first  order,  is  synchronous 
with  neither  if  we  take  8^  =1.  The  case  given  by  Eq.  (2.71) 


I 


i 


$ 

1 


i 


s 

j 

! 

>- 

i 


in  which  the  zero  and  first  orders  are  coupled,  meets  all  criteria. 


1*9 


It  is  satisfying  that  this  equation  is  equivalent  to  Eq.  (2.1*1*) 
derived  earlier  in  describing  Bragg  scattering. 

With  only  MQ^^(n)  and  M^P^(r))  nonvanishing,  Eq.  (2.1*6) 
provides  the  system  of  coupled  mode  equations 


2 


dn 


(2.72) 
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+  M  (p) 
o 


0 


(2.73) 


Solutions,  subject  to  the  boundary  condition  (2.1*8),  are 


Mq(p)  =  D  cos  (|  n)  (2.7»*) 

Mx(p)  =  — D  sin  (|  n)  (2.75) 


In  terms  of  the  original  quantities,  the  intensities  at  x  =  L  are 
given  by 


k 
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I {p) 
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■  rpf<j) 

°f 

n<p) 

L) 

(2.76) 
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h(p) 

-  1)2  s1"2  V'1’ 

°f 

n(p) 

L) 

(2.77) 

The  approximation  w  «  uj  was  used  in  writing  Eqs.  (2.76)  and 

o  n 

(2.77).  Again,  in  the  limiting  case  of  a  bulk  wave,  the  solutions 

agree  with  previously  published  results.  '  With  sin  (a)  2l  a  for 

3mall  a,  Eq.  (2.77)  reduces  in  the  case  of  low  acoustic  power  to 

k 


{p)  = 


V"  -  ft  Vy> 


wf 

ITpT 


l] 


(2.78) 
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Comparison  of  Eqs.  (2.78)  and  (2.6?)  shows  that  the  first  order 
diffraction  profile  is  independent  of  whether  the  diffraction 
occurs  in  the  Bragg  or  Raman-Nath  regime. 

As  was  previously  indicated,  at  high  acoustic  frequencies 
it  is  possible  that  the  polarization  of  the  Bragg  diffracted  light 
is  rotated  from  that  of  the  zero  order.  The  procedure  for  deter¬ 
mining  the  intensity  of  the  diffracted  order  is  exactly  as  that 
employed  above,  except  that  in  this  case  R  (y)  is  replaced  by 

R  (y).  Thus,  when  polarization  rotation  occurs, 

PQ 


Io(r>  =  D* 


n 

k 


(p) 


L) 


T  ( Q. )  «*2  .  2  / 1  -j  /  \  _ f  T  \ 

h  =  D  sin  (¥  Rpq(y)  L) 


(2.79) 


(2.8o) 


2.6  Numerical  Results 

The  theory  of  the  acousto-optic  interaction  in  surface 
acoustic  waves  has  been  developed  in  the  preceding  sections  of 
this  chapter.  In  this  section  we  evaluate  the  obtained  expres¬ 
sions  for  the  special  case  of  a  z-propagating  Rayleigh  wave  on 

y-cut  lithium  niobate.  The  form  of  the  Rayleigh  wave  as  given 

( 68 ) 

by  Spaight  and  Koerber 


is 
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-J°I  V  -Jot2  V  -Ja3kay,  **J(“at-V) 

u2  =  B&  (e  +  m21  e  +  m^  e  )  e 

„  -IV,'  -Jt,2  “ay  ^  ■'•jV,  J(W> 

u3  =  Ba  N^e  n21  6  +  n31  e  )  e 


-ja^y  -Jot2kay  -ja,koy  j(wot-k„z) 


_  _  ,  n,  l  a-'  o  2  a  A  3  a'  .  - '  a  "a 

e2  ~  Ba  ka  °(e  +  a21  6  +  a3l  e  5  6 


R  v  p,  "Ja2kay  -Ja3kay,  J(V~kaz) 

e3  =  Ba  ka  P(e  +  b21  e  +  b3l  e  }  e 


el  =  U1  =  °*  y  -  0 


(2.81) 


Parameters  given  in  these  equations  are  provided  in  Table  2.1. 

Table  2.2  below  provides  numberical  values  for  the  para¬ 
meters  which  optically  characterize  lithium  niobate.  A  few  comments 
are  in  order  regarding  several  of  the  constants  listed.  First,  the 

usual  convention  has  been  employed  in  which  certain  sets  of  sub- 

(l6 ) 

scripted  indices  are  given  in  a  contracted  notation.  The  sets 

i,j  and  k£  appearing  in  PiJkJi  and  riJm  (i,j,k,  =  1,2,3)  are 

replaced  by  a  single  index  with  values  1,2, 3,..., 6.  Thus  P  ^ 

becomes  PT  T  and  r . .  becomes  rT  . 

IJ  ijm  Im 


i,j  -  11  1  =  1 

22  2 

33  3 

23  or  32  k 

13  or  31  5 

12  or  21  6 


52 


Symbol 


a,. 


a„ 


0 

P 


m. 


21 


m 


31 


‘21 


31 


21 


31 


21 


31 


TABLE  2.1 


Rayleigh  Wave  Parameters  for  YZ  LiNbO^^ 


Value 


3.5105  x  10 J  m/sec 
0.06462  -  3  0.1217 
0.3804  -  j  1.0379 
-0.3960  -  J  0.7738 


-0.05180  +  j  0.09717 


(-0.1826  +  j  0.1402)  x  109 

(-1.5202  -  J  0.6931)  x  109 

-0.02637  +  3  .3719 

0.02740  -  3  0.5141 


-1.9258  +  j  5-0939 
-4.U64  -  3  1.4533 
-32.435  -  3  22.191 


27-233  +  j  19-453 
-3.6274  -  j  3.2912 
-0.05632  +  3  5-3042 


Description 
Phase  Velocity 


Decay  Constants 


Wave  Amplitude 
Coefficients 


B  .j."—  '  A*'i  >.  ' 


TABLE  2.2 


Optical  Properties  of  LilfbO^ 


Symbol 


Value 


Description 


e°  = 

£° 

5.24 

11 

22 

c° 

4.84 

33 

P  .  = 

P„„ 

-.02 

11 

22 

P  „  = 

P„ 

.08 

12 

21 

P„  „  = 

P„„ 

.13  e' 

13 

23 

Pl4  = 

~P24  =  P65 

-.08 

p„  = 

P 

.17 

31 

32 

P__ 

.07 

33 

P4l  = 

_P42  =  P56 

-.15 

P44  = 

P55 

.12 

33 


r22  rl2  r6l 


r13  =  r23 


r51  "  rU2 


J77  * 


Unperturbed  Relative 
Permittivity ^ ) 


Elasto-Optic  Tensor 
Element  ^ ' 


30.8  x  10~2  ej(‘*125lT)  m/V*  Electro-Optic 

Tensor  Element'1-5' 

3.4  x  10~12  m/V 
8.6  x  10"12  m/V 
28  x  10~12  m/V 


modified  value 


The  values  employed  for  the  elasto-optic  coefficients  at  constant 

electric  field  and  electro-optic  coefficients  at  constant  strain 

(27)  (73) 

were  those  given  by  Pinnow  and  Kaminow  respectively.  It 

was  found  in  calculating  Ae.  from  these  constants  using  Spaight 

and  Koerber  data  for  the  Rayleigh  wave,  numerical  discrepancies 

(19  7l*) 

existed  with  Ae^  calculated  using  Dixon's  *  effective  elasto- 
optic  coefficients  which  include  both  electro-  and  elasto-optic 
contributions.  This  discrepancy  was  minimized  by  introducing  a 
phase  factor  on  some  of  the  constants,  which  is  reflected  in 
Table  2.2.  This  resulted  in  good  agreement  with  experimental 
results  of  Kramer. 


With  reference  to  Tables  2.1  and  2.2  and  Eqs.  (2.10) 


and  (2.38), 


(2.82} 


assuming  =  6328  A.  This  frequency  was  not  approached  in  any 
of  the  experiments  reported  here. 

Table  2.3  gives  the  calculated  values  of  the  coefficients 
appearing  in  Eq.  (2.18).  Coefficients  for  the  elasto-optic, 
electro-optic,  and  combined  effects  are  provided.  Based  on  these 
coefficients,  the  function  (y) ]  is  calculated.  The  intensity 
profile  for  light  in  the  first  diffraction  order  is  proportional 
to  this  function,  assuming  low  acoustic  power.  Results  are  shown 


graphically  in  Figs.  2.3,  2.1*,  and  2.5. 
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Figures  2.3  and  2.k  have  obvious  implications  to  inte¬ 
grated  optics  structures  fabricated  on  lithium  niobate.  Those 
devices  employing  z-polarized  light  require  very  shallow  waveguides, 
less  than  oneeighth  of  an  acoustic  wavelength,  for  example.  In 
this  category  are  all  those  devices  fabricated  by  out-diffusion 
of  LigO  from  the  surface,  since  this  process  chiefly  effects  the 
extraordinary  refractive  index  and  therefore  provides  waveguiding 
for  z-polarized  light.  On  the  other  hand,  for  devices  which  can 
make  use  of  y-polarized  light  a  deeper  guide,  extending  to  about 
one  third  of  an  acoustic  wavelength,  would  take  better  advantage 
of  the  available  acoustic  power.  A  buried  waveguide  would  be 
most  advantageous  in  this  case. 

Figure  2.6  shows  the  integrating  effect  of  a  finite  beam 
width,  as  discussed  in  connection  with  Eq.  (2.68).  The  integral 

oo 

5  rect  (h-  I  1  (L,y)dy  ,  divided  by  (D  B  k  )2  , 

_oo  v  n  a  a 

is  shown  in  the  figure  for  I  =  2,3  (contracted  notation).  The 
solid  line  indicates  z-polarized  light;  the  dashed  line  is  for  the 
y-polarization.  The  trapezoidal  integration  technique  was  used  in 
generating  these  curves.  Sixteen  intervals  were  utilized  spanning 
a  total  width  w  =  .255  A«  As  noted  earlier,  the  diffraction  effect 
is  more  pronounced  near  the  surface  for  z-polarized  light.  Further 
from  the  surface  the  diffraction  of  y-polarized  light  is  stronger. 
These  calculated  results  are  superimposed  upon  the  experimental 
results  of  Kramer  in  Figs.  2.  7  and  2.8,  with  good  agreement.  In 


6l 


Fig.  2.7  it  was  necessary  to  translate  Kramer's  curve  approxi¬ 
mately  .08  A  toward  the  surface  in  order  to  have  the  maxima  of 
the  two  curves  coinciding.  This  is  probably  within  his  experi¬ 
mental  resolution  in  determining  the  crystal  surface  position 
with  the  optical  probe  he  used. 


CHAPTER  3 


ACOUSTO-OPTIC  INTERACTION  IN  SURFACE  ACOUSTIC  WAVES: 

DEVICE  APPLICATIONS 

3. 1  Unified  Description  of  Acousto-Optic  Devices 

In  this  section  we  develop  a  description  of  acousto¬ 
optic  devices  which  can  he  applied  to  the  operational  analysis 
of  such  varied  devices  as  modulators,  signal  processors,  and 
image  scanners.  Ihe  description,  which  is  appropriate  to  both 
bulk  and  surface  wave  devices,  is  developed  for  an  interaction 
in  the  Raman-Nath  regime. 

Consider  the  one- dimensional  optical  system  shown  in 
Fig.  3.1.  A  uniform  plane  wave  is  amplitude  modulated  by  the 
transmittance  function  T(^)  which  is  imaged  via  lens  into 
the  z-plane,  where  it  appears  as  T(z).  Depending  on  the  place¬ 
ment  of  L1  the  imaging  process  could  involve  magnification. 

Since  T(z)  is  spatially  limited  we  can  consider  an  aperture  of 
width  &  £  d,  positioned  in  the  z-plane  of  the  diffraction 
limited  system.  Voltages  f^(t)  and  f^ (t )  amplitude  modulate  rf 
carriers  at  the  acoustic  frequency  us  which  excite  transducers 
at  opposite  ends  of  the  ultrasonic  delay  line.  It  is  assumed 
throughout  that  the  time  rates  of  change  of  f^(t)  and  f2(t)  are 
small  compared  to  that  of  the  rf  carrier.  These  signals  produce 
traveling  optical  phase  gratings  which  diffract  the  incident 
light.  The  Fraunhofer  diffraction  pattern  is  formed  in  the  z'- 
plane,  which  is  the  focal  plane  of  lens  L^. 
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The  amplitude  distribution  of  the  diffracted  light  is 


given  by  the  diffraction  integral 


.  *  f  £/2 

A(6)  =  CeJa,t]  T(z) 


J  [y-^-  sin  0  +  V1  (v&t  - 


z  -  — )  sin  (oj  t  -  k  z ) 
2  a  a 


+  V„  (v  t  +  z  +  ^)  sin  (u*  t  +  k  z)] 

2  a  2  a  a 

dz  ,  (3.1) 

in  which  C  is  a  normalizing  constant,  u)  is  the  frequency  of  the 
incident  light,  A^  is  the  free  space  wavelength  of  the  incident 
light,  the  angle  0  is  measured  from  the  x-axis,  v  is  the  acoustic 
velocity  given  by  oi  /k  ,  and  the  spatial  functions  V  and  V  are 

&  &  X  CL 

proportional  to  the  applied  voltages  f^(t)  and  f2(t).  In  particular, 
with  respect  to  a  common  point  of  insertion, 

VV>  -  'i<t>  (3.2) 

°f 

and 

V2(vat)  =  Z-UX  f  (t)  ,  (3.3) 

°f 

in  which  b  is  a  constant  which  takes  into  account  material  photo¬ 
elastic  constants  and  transducer  electro-mechanical  coupling. 

Although  Eq.  (3.1)  is  the  diffraction  integral  of  scalar  diffraction 
theory,  we  interpret  A(0)  as  a  component  of  the  electric  field  of 


the  diffracted  light. 
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Using  the  identity 

QO 

Rasing  =  Z  j  (a)  ejr6  ^ 


(3.M 


r  =  -°° 

in  which  J  is  the  order  Bessel  function  of  the  first  kind, 
r 


Eq.  (3.1)  becomes 


.  .  r  A/2 

A(0)  =  C  eJ“M  T(z)  %  S  J  [Vn(v  t-z-d/2)]  J  [V_(v  t+z+d/2)] 
-A/2  s  r  r  1  a  s  ^  a 


sin  0  z  +  (r  +  s)  w  t  +  (s  -  r)  k  z] 
A  a  a 


x  e 


dz 


(3.5) 


Assuming  uniform  convergence  of  both  series,  the  order  of  integra¬ 
tion  and  summation  can  be  exchanged.  Letting  n  =  r  -  s,  Eq.  (3.5) 
becomes 


A(6)  =  C  eJ<llt 


r  A/ 2 

n  T(z)  JjV,(v  t-z-d/2)]  J 

^  0  /o  r  1  a  r 

n  r  -A/2 


[V  (v  t  +  z  +  d/2)] 
n  2  a 


J[(2r-n)  to  t  + 

a 


x  e 


(jp^-  sin  0 
°f 


n  k  )z ] 
a 

dz 


(3.6) 


For  ideal  diffraction,  A  -*•  and  the  integral  has  nonvanishing 
values  only  for 


sin  0  -  n  k 

a 


0 


(3.7) 


which  means  that  light  occurs  only  at  discrete  angles  given  by 
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n  A 

°f 

sin  6  =  t -  .  (3.8) 

n  A 

In  this  case  the  light  in  the  z' -plane  appears  as  6-functions  at 
the  points 

n  F  A 

of 

z’  i  F  sin  6  =  - — -  (3.9) 

n  n  A 

In  the  physical  case  in  which  the  aperture  &  is  finite,  diffraction 
spreading  of  these  spots  occurs,  the  width  of  each  spot  given 
approximately  by  F  Aq  /£. 

The  amplitude  of  the  light  in  the  n^  diffracted  order 
is  given  by 

j[(2r-n)o)  +oi]t  *A/2 

A  =  CZe  a  \  T(z)  J_[V  (v  t- z- d/2] 

n  r  -1/2  1  a 


J  n  [V_(v  t+z+d/2)]  dz 
r  -  n  d  a 


(3.10) 


The  time  dependent  intensity  in  the  n 


th 


order  is  given  by 


I 

n 


A  A  * 
n  n 


J2(r-p)o)  t  Cl/ 2 


8  C‘  H  e  “  J  T(z)  J  [V  (v  t-z-d/2)]J  [V9(v  t+z+d/2 ) ]  dz 

r  p  -A/2  r  1  a  r-n  2  a 


pA/ 2 

x  \  T(z)  J  [V  (v  t-z-d/2)]  J  [V_(v  t+z+d/2) ]dz 
-A/2  ?  1  a  p-n  2  a 


(3.11) 


The  time-averaged  light  intensity,  taken  over  a  time  interval 
greater  than  several  rf  cycles  but  small  compared  to  the  time  rate 
of  change  of  f^(t)  and  fg(t),  is 
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P  r«./2  2 

I  =  E  I J  T(z)  J  [V  ( v  t-z-d/2 ) ]  J  [V- ( v  t+z+d/2 ) ]dz | 

n  r  =  _oc  -1/2  1  r-n  ^  a 


(3.12) 


To  relate  Eq.  (3.12)  to  the  results  of  Chapter  2,  we 


take  Vg  =  0,  and  T(z)  =  1,  with  the  result  that  Eq.  (3.12)  becomes 


I  =  C 
n 


2 


|\  J  [V  (v  t-z-d/2)]  dz 
-SL/2  n  1  a 


(3.13) 


In<i  ^total  Ll 


dy 


If  f^(t)  is  constant  with  time,  this  can  be  further  reduced  to 

I  =  C2  £2  J  2  (V  ) 
n  n  1 

th 

The  total  time-averaged  light  intensity  in  the  n  order  is 
obtained  from  the  results  of  Chapter  2  by  taking  the  integral 
over  y  of  Eq.  (2.65). 

'h  +  w/2 
w/2 

Equation  (3.15)  represents  the  total  light  in  the  nth  order 

at  x  =  L.  Parseval's  theorem  provides  that  this  is  the  same 

til 

intensity  that  would  appear  in  the  n  order  in  the  z' -plane. 

The  mean  value  theorem  for  integrals  guarantees  that  there  exists 

some  value  y  =  a  in  the  interval  h-w/2£y£h  +  w/2  such  that 

the  integral  of  Eq.  (3-15)  can  be  written 

k 

r  (i) 


(3.1M 


(3.15) 


total 


■  ',I>2  J„2  >2  "ii(a)  737  Ll 


(3.16) 


A  comparison  of  Equations  (3.11*)  and  (3.16)  provides 
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and 


(3.17) 


(3.18) 


From  Eq.  (3.18)  it  is  apparent  that  the  constant  b  appearing  in 
Eq.  (3-2)  is  given  by 


b 


1  Rii(ct) 

2 


(3.19) 


which  is  a  measure  of  the  effectiveness  of  an  applied  voltage 
signal  to  produce  a  perturbation  of  the  material's  permittivity. 
We  are  now  in  a  position  to  indicate  several  device  applications 
of  the  acousto-optic  interaction. 


Modulator 

The  modulator  is  the  simplest  acousto-optic  device.  By 
varying  the  input  voltage  f^(t),  the  time-averaged  light  intensity 
of  the  various  orders  is  controlled.  With  f2(t)  =  0  and  T(z)  =  1 

x  -u 

in  -  £/2  £  2  .1  the  intensity  in  the  n  order  is  given  by  Eq. 

(3-lU).  Using  Eq.  (3.2), 

In  =  iQ  Jn2  [K  fx(t)]  ,  (3.20) 

in  which  K  is  a  constant  given  by  2w  L  b/X  .  It  should  be  recalled 

°f 

that  Eq.  (3.lM  was  derived  under  the  assumption  that  was  constant 
in  the  optical  aperture  region.  Thus  we  require  that  the  aperture  £ 
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be  small  compared  to  the  spatial  rate  of  change  of  V^.  In  terms 
of  the  angular  frequency  of  the  highest  frequency  component  of 
fl(t),  a)  we  can  express  this  requirement  as 

n  v 

£  <  - -  ,  (3.21) 

<u  * 

max 

in  which  the  acoustic  wavelength  corresponding  to  w  is  at 

IucLX 

least  twice  £.  For  a  signal  with  w  =  2n  x  1  MHz,  this  requires 

nicix 

£  to  be  less  than  1.75  mm,  which  is  a  typical  beam  diameter  for  a 
low  power  HeNe  laser. 

Of  particular  interest  is  the  case  n  =  1  for  a  small 
signal  f^(t).  In  this  case, 

I  IQ  K2  [^(t)]2  .  (3.22) 

The  first  order  diffracted  light  is  proportional  to  the  power  of 
the  input  signal. 

Image  Scanner 

With  fg(t)  =  0,  the  configuration  shown  in  Fig.  3.1  can 

2 

be  used  to  transform  an  intensity  transmission  function  t(^)  =  |t(c)| 
into  a  voltage  for  the  purpose  of  facsimile  transmission.  From  Eq. 
(3.12), 

?  1/2  2 

I  =  CT  I  T(z )  J  [V,  (v  t  -  z  -  d/2)]  dz  .  (3.23) 

n  _l/2  n  1  a 

For  small  arguments  of  the  Bessel  function,  the  first  order  light 


is  given  by 


TO 


P  1/2  2 
I,  =  <T  |  T( z )  V.  ( v  t  -  z  -  d/2)  dz  |  .  (3.21*) 

1  -1/2  X  a 

Equation  (3.21*)  can  be  interpreted  as  the  square  of  the  convolution 

of  the  spatial  mask  T(z)  with  the  acoustic  signal  V^.  If  f^(t)  is 

a  very  short  rectangular  pulse  which  amplitude  modulates  the  rf 

carrier,  then  V  approximates  the  delta  function,  6(v  t-z-d/2). 

1  8, 

Thus  the  integral  of  Eq.  (3.23)  becomes 

\  =  C2  |T(vat  -  d/2) |2 

=  C2  t(v  t  -  d/2)  (3.25) 

As  a  function  of  time  the  intensity  replicates  the  spatial  infor¬ 
mation  t(^).  In  practice  need  not  be  made  small  in  amplitude. 

The  pulse  height  is  adjusted  to  maximize  the  first  order  Bessel 
function. 

Convolver 

It  is  apparent  from  Eq.  (3.11)  that  the  time  dependent 
intensity  in  each  diffraction  order  is  made  up  of  a  sum  of  temporal 
frequency  components  occurring  at  even  integral  multiples  of  the 
acoustic  rf  frequency.  Any  given  frequency  component  can  be  elec¬ 
tronically  isolated  by  proper  filtering  of  the  output  signal  of  a 
photodetector  placed  to  collect  the  light  in  a  single  diffraction 
order.  The  resultant  signal,  after  filtering,  provides  the  real¬ 
time  convolution  of  f^(t)  and  f2(t).  We  illustrate  this  for  the 
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case  of  light  in  the  zero  order  at  twice  the  acoustic  frequency, 
but  a  similar  analysis  can  be  applied  to  obtain  real-time  convolu¬ 
tion  in  any  order  at  any  even  multiple  of  acoustic  frequency. 

In  the  special  case  that  f^(t)  and  fg(t)  are  time 
limited  signals , 


fx(t),  f2(t)  =0  for  t  <  0,  t  >  d/va,  (3.26) 

and  T(z)  is  a  simple  window  function. 


f  1  1*1  <  d/2 

T(z)  =  J  (3.2?) 

^  0  i z |  >  d/2  , 

the  limits  of  integration  for  all  integrals  appearing  in  Eq.  (3.11) 
can  be  taken  from  - «°  to  00  without  affecting  the  value  of  the  inte¬ 
grals.  With  n  =  0  for  the  zero  order  diffracted  light,  and  (r-p)  =  1 
corresponding  to  the  component  at  twice  the  acoustic  frequency,  Eq. 
(3.11)  provides,  after  rearrangement  of  terms. 


(2<0  p  r°° 

IQ  =  2CT  {  J  J1[V1(vat-  z-  d/2) ]  J1[V2(v&t  +  z  +  d/2)]  dz} 


x(\  J  [V,(v  t- z-d/2)]  J  [V.(v  t  +  z  +  d/2 ) ]  dz 
v  _  ola  o  d  a 


Joo 

J2l'Vl(vat- z- d/2)]  J2[V2(vat  +  z  +  d/2)]  dz) 


+  (higher  order  terms) 


(3.28) 
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in  which  terms  involving  sixth  and  higher  order  products  have  teen 
neglected.  The  function  J-^V^)  J^(Vg)  ^Jo^Vl^  + 

J2(Vi>  j  (v2)]  is  plotted  in  Fig.  3.2  for  various  standing  wave 
ratios  V^/Vg.  For  small  arguments  of  the  Bessel  functions,  Eq. 
(3.28)  becomes 


<2»a) 


A  OO 

2t  |  C2  d  J  W-z-d/2)  V2(v&t  +  z  +  d/2)  dz  .  (3.29) 


With  the  substitution  v  T  =  v  t  -  z  -d/2,  Eq.  (3.29)  reduces  to 

a  a 

the  familiar  convolution  integral.  In  terms  of  the  voltage  sig¬ 
nals  f^(t)  and  f2(t). 


(2io  )  2  T2  2 

I  a  %  C2  d 

O  ~  y  2 


J ; 


)  fg(2t  -  x)  dT 


(3.30) 


An  interesting  feature  of  Eq.  (3.30)  is  the  time-compression 
factor  of  two  which  arises  from  the  fact  that  both  functions 
and  V2  are  propagating  in  opposite  directions  along  the  acoustic 
track.  In  the  usual  convolution  integral,  one  function  is  con¬ 
sidered  fixed  with  respect  to  the  coordinate  system  while  the 
other  is  shifted  with  respect  to  it. 

The  analog  convolver  is  a  powerful  signal  processor 
because  of  its  great  flexibility.  If  the  two  input  carrier  fre¬ 
quencies  are  not  identical,  and  if  f2(t)  is  an  even  function  of 
time,  then  the  modulating  factor  of  the  output  frequency 
2(u>  -  Aw  )  is  proportional  to  the  cross-ambiguity  function  of 

8l  81 

fj.  and  f2: 
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[2(w  -  Aw. )  ]  f*  "  jAw  T 

IQ  a  =j  f^T)  f2(r  -  2t)  e  a  dt 

W  — oo 


(3.31 


In  the  case  that  Aw  =  0,  Eq.  (3.31)  reduces  to  the  correlation 
a 

integral.  Equation  (3. 31)  also  provides  a  means  of  producing  the 
Fourier  Transform  of  the  function  f^(t),  by  taking  f2(t)  equal  to 
a  constant  and  sweeping  Aw  . 

Table  3.1  on  the  next  page  indicates  the  application 
capabilities  of  the  optical  system  of  Fig.  3.1.  In  the  table, 
a  6-function  optical  input  T(z)  refers  to  a  narrow  window  func¬ 
tion;  a  6-function  voltage  input  refers  to  a  short  rectangular 
pulse. 


3.2  Convolver  and  Related  Devices 
3-2.1  Introduction 

In  its  simplest  form,  a  communication  system  is  charac¬ 
terized  by  a  transmitter  end  where  information  is  encoded  and 
imposed  upon  some  carrier,  a  transmission  medium,  and  a  receiver 
end  where  the  information  is  abstracted  from  the  received  signal. 

A  major  consideration  in  all  real  communication  systems  is  the 
signal  to  noise  ratio  at  the  output  of  the  receiver.  From  a 
practical  point  of  view,  signal/noise  considerations  determine 
the  number  of  repeaters  required  along  the  transmission  path, 

and  thus  have  a  direct  bearing  on  the  cost  of  the  system. 

(75 ) 

It  is  well  known  that  the  signal  to  noise  ratio  with 
respect  to  power  at  the  output  of  a  receiver  is  given  by 
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S 

N 


t*  00 

1A  IT2  |  )  S{(d)  H(oj)  eJ  1  du>| 


f "  2 

N£^2tt  J  |H(a))|  du> 


(3.32) 


where  o>  represents  angular  frequency,  t  is  time,  S (oj )  is  the  fre¬ 
quency  spectrum  of  the  received  signal,  H(w)  is  the  trans¬ 
fer  function  of  the  receiver  which  is  assumed  linear  and  time 
invariant,  and  NQis  the  noise  spectrum  power  density,  which  for 
simplicity  has  been  assumed  flat.  For  a  given  signal,  this  ratio 
can  be  maximized  by  taking  the  receiver  frequency  characteristic 
to  be  of  the  form 


H(w)  =  aS*(u>)  e'ja)T  (3.33) 

* 

where  S  (w)  is  the  complex  conjugate  of  the  signal  spectrum,  and 
a  and  t  are  arbitrary  real  constants.  This  approach  to  matched 
filtering  has  of  course,  found  widespread  use  in  electrical  and 
optical  processing.  The  output  of  such  a  matched  receiver  to  a 
signal  of  frequency  spectrum  S(w)  is  given  by 


V 

o 


1/2tt 


p  00 

J  S(w)  S  (w)  e*"^aVt  aw 


(3-3**) 


where  for  simplicity,  we  have  taken,  a  =  1,  x  =  0.  It  is  easily 
shown  that  this  output  can  be  equivalently  written  as  an  auto¬ 
correlation  function  of  the  signal  s(t)  itself,  where 


s(t) 


l/2ir 


J  S(w)  e^wt  du 


(3.35) 
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P  CO 

S(w)  =  J  s { t )  e"^  dt  (3-36) 

.CD 

That  is 

A  co 

V  =  J  s ( t * )  s(t'  -  t)  dt  (3.37) 

—  00 

For  the  receiver  output  in  a  real  communication  system,  s(t')  in 
Eq.  (3.37)  would  be  replaced  by  s(t')  +  n(t'),  where  n(t’)  repre¬ 
sents  the  noise.  It  is  apparent  from  Eqs.  (3-3^)  and  (3.37)  that 
best  reception  requires  that  either  the  receiver  be  matched  in  its 
frequency  response  to  the  incoming  signal,  or  that  the  receiver  be 
capable  of  signal  correlation.  The  former  method  requires  that 
the  receiver  have  a  different  frequency  response  for  each  signal 
anticipated;  the  latter  method  requires  a  means  for  real-time 
correlation.  These  mathematically  equivalent  methods  have  very 
different  realizations.  The  experimental  demonstration  of  real¬ 
time  correlation  via  the  surface  acousto-optic  interaction  is  pre¬ 
sented  in  Section  3-2.2. 

If  s(t)  represents  a  transmitted  radar  signal,  then  the 
matched  filter  output  of  the  receiver  for  the  return  signal 
reflected  from  a  moving  target  is  given  by 

Vo  =  J  s(t')  s(t*  -  t)  ejAwt'  dt  .  (3-38) 

.00 

In  Eq.  (3.38),  which  represents  the  ambiguity  function  of  s(t),  the 
time  delay  t  is  related  to  the  target's  range,  and  the  doppler  shift 
Aw  is  related  to  the  time  rate  of  change  of  the  range. 
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The  use  of  waveforms  which  can  be  processed  at  the 
receiver  to  enhance  their  detectability  is  of  particular  impor¬ 
tance  in  modern  communications.  Each  data  element  in  such  a 
system  is  generally  encoded  into  a  symbol  from  sin  a  priori 
established  alphabet.  It  is  usual  that  the  transmitted  waveforms 
are  of  more  complex  form  than  the  baseband  data,  but  correlation 
can  be  performed  at  the  receiver  to  increase  detection  probability 
since  each  symbol  is  one  of  a  set  of  predetermined  waveforms.  The 
increased  detectability  is  achieved  at  the  expense  of  increased 
bandwidth  requirements. 

One  code  which  has  found  widespread  use,  especially  in 
multiple-access  time  slot  systems,  is  the  Barker  code.  This  code 
is  especially  attractive  since  the  autocorrelation  of  an  N-bit 
Barker  code  has  a  correlation  peak  which  is  N  times  the  magnitude 
of  the  side  lobes.  Barker  sequences  {q},  defined  through  the 
relation 


N-k-1  (  0 ,  N  -  k  even  with  k  #  0 

Z  qi  qi  +  k  =  1 
i  =  0  V  +  1 ,  N  -  k  odd 


(3.39) 


have  been  found  empirically  for  N  =  2,  3,  5,  7,  11,  13.  The 

individual  elements  may  take  on  values  +1  or  -1.  For  N  >  13,  no 
Barker  codes  exist.  It  is  to  be  mentioned  that  other  important 
codes  do  exist  and  have  been  of  use.  Of  special  interest  are  the 
pseudo- random  noise  code  and  linear  fm  chirp. 
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3.2.2  Experiment 

The  basic  optical  processor  is  shown  schematically  in 
Fig.  3-3.  Light  from  a  3  mW  HeNe  laser  is  brought  to  a  line  focus 
by  lenses  to  L^,  producing  a  n-  15  pm  wide  line  in  the  surface 
of  the  LiNbO^  delay  line.  Lens  L^  recollimates  the  light  and  Lg 
produces  the  Fraunhofer  diffraction  pattern  of  the  light  amplitude 
distribution  in  its  focal  plane.  Amplitude  modulated  signals  are 
applied  to  inputs  1  and  2,  and  the  final  optical  signal  is 
detected  by  a  photodiode  positioned  in  the  focal  plane  of  Lg 
where  it  collects  the  light  of  a  single  diffracted  order. 

The  output  of  the  photodiode,  after  electronic 
processing,  is  displayed  on  an  oscilloscope.  A  typical  processing 
scheme  is  shown  in  Fig.  3-*+.  Heterodyning  of  the  photodiode  sig¬ 
nal  and  subsequent  filtering  with  a  tuned  amplifier  isolates  the 
voltage  component  at  the  desired  frequency. 

Y-cut  Z-propagat ing  LiNbO^  was  employed  for  the  acousti¬ 
cal  delay  line  with  the  interdigital  transducers  matched  to  50  £2 
at  the  center  frequency.  Experiments  were  performed  at  several 
frequencies,  results  being  given  for  ^5  MHz.  The  ^5  MHz  delay 
line  had  a  bandwidth  of  9  MHz  and  an  insertion  loss  of  approxi¬ 
mately  6  db  per  transducer.  The  2  cm  separation  of  the  trans¬ 
ducers  corresponded  to  a  total  delay  of  5-8  ysec.  The  transducer 
apertures  were  2  mm  wide.  Diffraction  was  in  the  Raman-Nath 
regime.  Convolution  was  obtained  in  zero,  first,  second,  third, 
and  fourth  diffraction  orders. 


I'll!.  >•')  I’l'Vica)  H-occuaor  lor  K.nl-Tl.u  Convolution 


OSCILLOSCOPE 
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Fig.  3.4  Electronic  Processing  for  Convolution.  Electronic 
heterodyning  isolates  the  component  of  the  photodiode 
signal  at  twice  ths  acoustic  frequency. 
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Figure  3.5  indicates  results  of  the  different  methods 
of  obtaining  the  output  of  the  optical  convolver.  There  is  a 
choice  as  to  which  diffracted  order  in  which  to  position  the 
photodiode,  as  well  as  a  choice  of  how  to  process  the  signal  from 
the  photodiode.  The  input  square  pulses  of  500  mW  power  and 
2  ysec  duration  which  were  applied  to  both  transducers  are  shown 
in  trace  3.5(a).  In  3 - 5 ( ) >  on  the  same  time  scale,  is  the  total 
intensity  detected  in  the  zero  order,  displayed  on  a  high  fre¬ 
quency  oscilloscope  after  amplification.  The  signal  is  composed 
of  a  dc  component  and  also  high  frequency  components  containing 
the  desired  autoconvolution.  The  component  at  2u>a  is  obtained 
by  heterodyning  and  filtering  and  is  displayed  in  3.5(c),  which 
is  the  I.F.  signal  of  the  tuned  amplifier.  Figure  3.5(d)  shows 
the  tuned  amplifier  video  output,  which  is  the  desired  envelope 
of  the  autoconvolution  signal  at  2 u)  in  the  first  order. 

It  should  be  noted  in  these  pictures  that  following 
the  detected  signal  in  time  is  a  smaller  signal  which  is  unwanted. 
This  signal  is  due  to  the  reflection  of  the  acoustic  wave  upon 
encountering  the  interdigital  transducer  at  the  opposite  end  of 
the  delay  line.  A  method  was  developed  which  completely  eliminated 
this  reflection  problem. 

In  this  method  the  surface  acoustic  waves  do  not  super¬ 
impose,  but  run  parallel  to  each  other  by  offsetting  one  trans¬ 
ducer  slightly  from  the  other  as  shown  in  Fig.  3.6.  For  bulk  waves 

( 76 ) 

this  configuration  has  been  treated  theoretically  '  as  well  as 
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experimentally. 
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Fig*  3*5  Experimental  Results  of  Real-Time  Convolution 
(a)  Rectangular  pulse  which  is  applied  to  both  inputs 
(15Vp-p) ,  (b)  Total  intensity  detected  in  the  zero  order, 
(c)  Autoconvolution  signal  at  twice  the  acoustic  frequency 
detected  in  zero  order,  (d)  Envelope  of  autoconvolution 
at  twice  the  acoustic  frequency  in  the  first  order,  and 
(e)  Envelope  of  auto  convolution  at  four  times  the  acoustic 
frequency  in  the  first  order.  Horizontal  scale=  2p  sec/div. 
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(a) 


Pig.  3.6  Overlapping  and  Parallel  Aeouatic  Track* 

(a)  In  usual  acoustic  delayline  configuration  acoustic 
tracks  overlap.  Reflections  occur  at  opposite  transducer. 

(b)  In  parallel  acoustic  track  configuration  both  tracks 
can  be  terminated  with  black  wax  to  eliminate  reflections. 
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Pig.  3.9  Experimental  Arrangement  for  Convolver  Data 


Figure  3.7  shows  results  vising  the  parallel  track  con¬ 
figuration.  For  this  LiNbO^  sample  the  delay  line  length  was 
less  than  the  previous  2  cm,  being  now  1.5  cm.  The  absence  of 
the  autoconvolution  of  the  reflected  wave  is  readily  apparent. 
Trace  3.7(a)  is  the  rectangular  input  pulse,  3.7(b)  the  total 
intensity  in  the  zero  order,  3.7(c)  and  (d)  the  autoconvolution 
and  envelope  at  2a>  in  the  zero  order,  3.7(e)  and  (f)  the  auto- 
convolution  at  2u )  in  the  first  order,  and  the  envelope  of  the 
autoconvolution  at  2 u  in  the  first  order,  respectively. 

Figure  3.8  shows  results  for  the  autoconvolution  of  a 
signal  composed  of  two  rectangular  pulses.  Shown  in  the  figure 
are  (a)  the  input  signal,  (b)  the  autoconvolution  detected  at 
2(i)  in  the  zero  order,  and  (c)  the  autoconvolution  envelope. 

For  a  rectangular  acoustic  pulse  of  fixed  length,  it 
is  expected  that  the  peak  height  of  the  triangular  autoconvolu¬ 
tion  depends  on  the  magnitude  of  the  acoustic  signal.  For  an 
ideal  convolver,  the  peak  height  would  be  proportional  to  the 
input  signal  voltage.  However,  as  was  shown  in  Section  3.1, 
this  situation  only  pertains  in  the  small  signal  limit  for  the 
acousto-optic  convolver.  For  larger  input  signals  the  convolu¬ 
tion  height  follows  the  curve  of  Fig.  3.2.  From  Eq.  (3.18)  it 
is  also  apparent  that  the  convolution  height  is  sensitive  to 
polarization. 

The  experimental  arrangement  of  Fig.  3.9  was  used  to 
measure  the  peak  convolution  as  a  function  of  the  signal  power 


Fig*  3*7  Experimental  Results  of  Real-Time  Convolution 
Using  Parallel  Tracks. 

(a)  Rectangular  pulse  which  is  applied  to  both  inputs 
(15Vp-p)(b)  Total  intensity  in  the  zero  order,  (c)  Auto¬ 
convolution  and  (d)  Envelope  in  the  zero  order  detected 
at  twice  the  acoustic  frequency  (e)Autoconvolution  and 
(f)  Envelope  at  twice  the  acoustic  frequency  in  the  first 
order.  Horizontal  Scale=  2y  sec/div. 


L*. 


85 


Fig.  3*8  Experimental  Results  of  Real-Time  Convolution 
Using  Parallel  Tracks. 

(a)  Input  signal,  (b)  Autoconvolution  detected  at  twice 

the  acoustic  frequency  in  the  zero  order,  (c)  Autoconvolution 

envelope. 


PEAK 
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of  the  inputs.  The  acousto-optical  convolver  was  exactly  as  that 
shown  in  Fig.  3.3,  except  that  a  polarizer  was  placed  before  the 
LiNbO^  crystal  and  an  analyzer  after.  This  was  necessary  in  order 
to  take  separate  measurements  for  y-  and  z-polarization  since  the 
HeNe  laser  which  was  used  provided  randomly  polarized  light.  To 
increase  the  light  throughput,  the  analyzer  was  removed  after 
verifying  that  polarization  rotation  did  not  occur.  The  acoustic 
input  signals  were  single  rectangular  pulses  in  which  the  pulse 
length  was  not  varied.  Data  was  taken  in  different  diffraction 
orders,  as  well  as  for  different  frequency  components  of  the 
photodiode  voltage.  All  data  shown  has  been  compensated  for  the 
slight  nonlinearity  of  the  detection  system,  including  photodiode, 
mixer,  and  tuned  amplifier. 

Figure  3.10  shows  peak  convolution  voltage  as  a  function 
of  input  power  for  identical  rectangular  inputs.  The  graph  is  of 
the  first  order  diffracted  light  for  the  convolution  signal  being 
detected  at  twice  the  acoustic  frequency.  The  vertical  scale  of 
this  graph,  as  with  the  following  graphs,  has  been  normalized  to 
the  maximum  signal;  it  is  the  shape  of  the  curve  that  is  of 
interest.  The  horizontal  scale  is  of  input  power  on  a  decibel 
scale  in  which  0  db  corresponds  to  a  lw  p-p  input  divided  equally 
between  the  two  50  ^-matched  interdigital  transducers. 

For  both  z-polarized  and  y-polarized  light  the  convolu¬ 
tion  signal  increases  with  input  power,  attains  a  maximum,  and 
then  decreases  with  increasing  acoustic  power,  although  the  maxima 


occur  at  different  acoustic  povers  for  the  tvo  polarizations. 
Because  one  would  want  to  operate  the  device  on  a  monotonically 
increasing  portion  of  the  curve,  the  location  of  the  maximum  is 
important . 

Figure  3.11  is  again  for  the  first  order  diffracted 
light,  in  which  the  convolution  signal  appearing  in  the  photo¬ 
diode  voltage  at  four  times  the  input  acoustic  frequency  is  now 
detected.  As  in  the  previous  case  the  y-polarized  convolution 
signal  attains  its  maximum  at  a  lower  power  level  than  the  z- 
polarized  one,  all  other  factors  being  equal.  The  families  of 
curves  shown  in  Figs.  3.12  and  3*13  were  generated  by  taking 
the  inputs  to  be  of  unequal  magnitude.  The  peak  convolution 
signal  is  plotted  as  a  function  of  the  magnitude  of  the  larger 
input.  The  top  curve  in  each  graph  is  for  the  case  in  which  one 
input  is  3  db  down  from  the  other  input.  The  second  curve  is  for 
6  db  down,  and  the  bottom  curve  for  9  db  down.  Both  Figs.  3.12 
and  3-13  are  for  the  first  order  diffracted  light  and  for  the 
convolution  detected  at  twice  the  acoustic  frequency.  In  both 
graphs  the  maxima  for  the  three  cases  occur  at  the  same  point, 
although  this  point  depends  on  polarization. 

Figure  3.1*+  is  for  the  convolution  detected  at  twice 
the  acoustic  frequency  in  the  zero  order  light.  The  inputs  were 
equal  amplitude  rectangular  pulses.  Again  the  y-polarized  signal 
attains  its  maximum  before  the  z-polarized. 


PHASE-LOCKED 
SECOND  HARMONIC 
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Figure  3.15  shows  the  input  signal  and  detected  auto¬ 
correlation  for  a  single  rf  pulse,  a  3  tit  Barker  code,  7  tit 
Barker  code,  and  13  tit  Barker  code.  The  disproportionate  side- 
lote  level  in  the  last  case  was  due  to  an  electronic  difficulty- 
in  the  generation  of  the  code.  In  obtaining  these  results  a  100 
MHz  acousto-optic  convolver  was  employed,  using  the  side-entry 

interaction  configuration.  The  light  source  was  an  argon  laser 
o 

operating  at  51^5  A.  The  Barker  code  was  impressed  upon  the  rf 
carrier  by  introducing  a  phase  shift  of  it  in  the  electrical  cw 
signal  used  to  launch  the  acoustic  wave.  In  this  way  +1  was 
distinguished  from  -1.  In  detecting  the  autocorrelation  of  the 
Barker  code,  the  electronic  arrangement  of  Fig.  3.16  was  employed, 
in  which  the  second  harmonic  reference  signal  used  in  hetero¬ 
dyning  was  phase  locked  to  the  carrier  of  the  reference  code. 
Because  the  Barker  code  is  not  symmetric  with  respect  to  reversing 
the  order  of  the  bits,  it  was  necessary  to  send  the  time-inverted 
signal  from  the  correlating  end  in  order  to  perform  the  correlating 
function  with  the  convolver. 

The  autocorrelation  of  a  linearly  fm  chirped  signal  is 
shown  in  Fig.  3.17.  The  2  ysec  pulse  swept  at  a  rate  of 
5  MHz  sec  ^  provided  the  time  compression  factor  of  20  evidenced 
in  that  figure. 

As  was  pointed  out,  under  certain  conditions  the  analog 
convolver  can  be  used  to  generate  the  ambiguity  function  of  time 
limited  signals.  Figure  3.18  shows  the  electronic  arrangement  for 
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Fig  3*15  Autocorrelation  of  100MHz  Signals.  Result  for 
Barker  encoded  data. 

(a)  Input  rf  signal  and  (b)  Autocorrelation  of  3  bit  . 
Barker  code  (ly  sec/div.)  (c)  Autocorrelation  of  7  bit 
Barker  code.  (2 y  sec/div.)  (d)  Autocorrelation  of  13  bit 
Barker  code  (2p  sec/div.) 
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the  ambiguity  function  generation.  Figure  3.19  shows  results  for 
(a)  a  rectangular  pulse  and  (b)  a  signal  consisting  of  two  rectangu¬ 
lar  pulses.  The  multi-exposure  photographs  provide  a  three- 
dimensional  plot  of  the  ambiguity  function  as  a  function  of  t 
and  Aw. 

3.2.3  Discussion 

The  theoretical  treatment  of  the  acousto-optic  inter¬ 
action  of  the  previous  chapter  is  related  to  the  device  descrip¬ 
tion  employed  in  this  chapter  through  Eq.  (3.18)'.  It  is  of 
interest  to  relate  these  two  descriptions  of  the  interaction 
to  the  experimental  results  of  the  preceding  section. 

A  comparison  of  Figs.  3.2  and  3.11*  provides  an  experi¬ 
mental  determination  of  the  proportionality  constant  b  appearing 
in  Eqs.  (3.2)  and  (3.3).  From  Fig.  3.2,  for  the  case  in  which 
=  V  ,  the  maximum  convolution  signal  occurs  for  =  1.1.  In 
the  experimental  result.  Fig.  3.11*,  the  maximum  for  y-polarization 
corresponds  to  -U  db  from  the  f.l  V  point.  This  corresponds 

O 

to  a  voltage  signal  of  5-8  V.  Equation  (3-2),  with  X  =  6328  A 

-6  -1  ^ 

and  L  =  2  mm,  provides  the  result  b  =  9.6  x  10  V-  .  Similarly, 
for  z-polarization,  b  =  7.8  x  10  ^  V  ^ 

An  experimental  determination  can  also  be  made  for 
the  effective  phase  grating  R^(ot)  at  this  acoustic  power 

level.  With  =  1.1,  Eq.  (3.18)  provides  R^ot)  =  2.5  x 

_U 

10  ,  for  example.  These  numerical  values  should  only  be 


93 


Fig.  3.1?  Autocorrelation  of  Linear  FM  Chirp. 

(a)  Frequency  swept  input  pulse  (2ysec/div.)  (b)  Auto¬ 
correlation  (2psec/div.)  (c)  Autocorrelation  on  expanded 
scale  ( .1^  sec/div. ) 


Electronic  Arrangement  for  Ambiguity  Function 
Generation 


Fig.  3*19  Experimental  Results  for  Ambiguity  Function. 

(a)  Multiexposure  photograph  showing  slices  of  the  ambiguity 
function  surface  of' a  3psec  pulse,  (b)  Ambiguity  function 
for  a  signal  comprised  of  two  rectangular  pulses  each 
0.8y  sec  and  separated  by  approximately  1.5/“*  sec. 
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considered  estimates,  since  the  Bessel  function  expansion  Eq. 

(3.11)  was  truncated  in  obtaining  Fig.  3.2. 

In  Figs.  3.10-  3.1^  the  signal  associated  with  the  y- 
polarized  light  obtained  its  maximum  value  at  a  lower  acoustic 
power  level  than  that  associated  with  the  z-polarized  light. 

With  reference  to  Fig.  2.6  we  can  infer  that  the  optical  beam 
was  focused  one  third  of  an  acoustic  wavelength  or  deeper 
beneath  the  crystal  surface,  since  the  interaction  is  stronger 
there  for  y-polarized  light  than  for  z-polarized. 

From  an  operational  point  of  view,  several  comments  are 
in  order  regarding  the  use  of  the  devices  demonstrated  for  actual 
signal  processing.  Two  important  measures  are  time-bandwidth 
product  and  dynamic  range. 

The  time-bandwidth  product  is  given  by  the  product  of 
transducer  bandwidth  with  acoustic  delay  line  transit  time.  The 
former  is  proportional  to  transducer  center  frequency  and  inversely 
proportional  to  the  number  of  transducer  finger  pairs.  The  U5 
MHz  device  reported  had  a  time  bandwidth  product  of  52;  the  100 
MHz  device  had  a  time  bandwidth  product  of  200. 

For  the  processing  of  sophisticated  signals  comprised 
of  lengthy  algebraic  codes,  a  large  time-bandwidth  produc4  is 
required.  The  limitations  imposed  in  achieving  a  large  bandwidth 
product  are  primarily  technological  and  not  fundamental.  High 
frequency  transducers  are  required  and  long  crystals  of  good 
quality.  With  presently  available  technology,  time-bandwidth 
products  of  several  thousand  could  be  achieved. 
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The  dynamic  range  of  the  acousto-optical  convolver  was 
United  in  the  experiment  by  the  noise  level  of  the  I.F.  amplifier 
and  mixer  combination  employed.  Any  signal  had  to  be  above  approxi¬ 
mately  1  pV.  The  current  from  the  photodiode  is  given  by  3H,  where 
H  is  the  light  irradiance  and  3  is  the  diode's  response.  For  the 
HP  U203  photodiode  used,  3  is  1.0  pA  per  mW/cm  of  irradiance. 
Assuming  a  maximum  irradiance  of  1  W/cra  in  any  practical  experi¬ 
ment,  the  resulting  current  through  a  50  load  would  result  in  a 
voltage  of  50  mV.  On  this  basis,  the  maximum  obtainable  dynamic 
range  in  any  practical  experiment  would  be  20  log  (50,000)  =  91*  dB. 
The  dynamic  range  evidenced  in  Figs.  3.10-3-1^  is  substantially 
less  than  this  estimate  due  to  the  fact  that  in  the  experimental 

arrangement  employed  in  obtaining  this  data,  the  irradiance  never 

2 

approached  1.0  pA  per  mW/cm  . 

While  the  detectability  of  a  signal  in  noise  is  enhanced 

by  using  a  matched  filter  receiver,  it  should  be  noted  that  any 

physical  realization  of  a  real-time  correlator  will  have  its  own 

noise  contributions.  In  the  case  of  the  acousto-optic  convolver/ 

correlator,  the  chief  noise  contributions  are  those  associated  with 

the  light  detection.  An  especially  attractive  feature  of  the 

acousto-optic  signal  processor  presented  is  that  it  inherently 

employs  optical  heterodyning,  which  permits  quantum  limited 

(lU) 

detection,  thereby  minimizing  the  processor  noise  contribu¬ 
tion.  The  requirements  placed  on  the  photodetector  are  that  it 
should  have  a  linear  response  and  that  it  must  have  good  frequency 


response.  It  must  "be  capable  of  detecting  signals  of  at  least 
twice  the  acoustic  frequency. 

A  particular  feature  of  the  technique  used  to  generate 
the  ambiguity  function  should  be  mentioned.  A  separate  convolu¬ 
tion  operation  must  be  performed  for  each  discrete  value  of  the 
delay-time  parameter.  Thus  actual  radar  system  implementation 
would  require  a  signal  repeater,  which  would  sacrifice  real-time 
capability,  or  an  array  of  parallel  processors. 

Of  practical  importance  applicable  to  the  optical  system 
alignment  of  all  the  devices  discussed  in  this  chapter  is  the  fact 
that  the  photodetector  need  not  be  placed  exactly  at  the  diffrac¬ 
tion  plane,  nor  need  it  have  a  flat  surface.  It  is  spatially  inte¬ 
grated  intensity  that  is  measured,  and  the  total  power  of  the 
light  is  the  same  whether  measured  in  the  diffraction  plane  or 
behind  it,  as  long  as  the  intensity  pattern  does  not  spread  wider 
than  the  active  area  of  the  detector. 

3. 3  Optical  Communication  System 
3-3.1  Introduction 

To  use  the  convolver/correlator  described  in  the  Section 
3.2  in  an  optical  communication  system,  it  is  necessary  to  receive 
the  transmitted  data  at  the  receiver  end,  convert  this  to  an  elec¬ 
trical  signal,  vise  this  signal  as  one  input  to  the  acousto-optic 
correlator,  and  supply  the  desired  reference  code  to  the  other 


Far  more  desireable  would  be  to  use  the  light  directly 


as  one  of  the  correlator  inputs.  The  possibility  of  such  an 
arrangement  is  suggested  by  the  parallel  track  configuration  shown 
in  Fig.  3-6,  in  which  the  acoustic  tracks  are  nonoverlapping. 

A  natural  extension  of  the  parallel  track  configuration 
would  be  to  separate  the  two  tracks  completely,  using  one  for 
light  modulation  at  the  transmitter  end  of  a  communication  system, 
and  using  the  second  acoustic  track  for  correlation  of  the 
received  light  signal  at  the  receiver  end.  It  is  a  requirement  of 
this  particular  scheme  that  the  transmission  medium  preserve  the 
phase  front  of  the  transmitted  light.  Such  a  system  is  applicable 
where  the  optical  transmission  line  can  be  modeled  as  a  series  of 
lenses,  as  indicated  in  Fig.  3.20(a).  Under  this  assumption,  Fig. 

3.20(b)  depicts  the  communication  system  implementation. 

3*3.2  Experiment 

The  arrangement  shown  in  Fig.  3.21  was  used  to  approxi¬ 
mate  an  actual  communication  system.  A  uniform  argon  laser  beam 
was  brought  to  a  line  focus  in  the  transmitter  acousto-optic 
modulator.  Lens  L^,  placed  intermediately  between  the  two  LiNbO^ 
acousto-optic  modulators  and  twice  its  focal  length  from  both  pro¬ 
duced  a  real  inverted  image  with  unit  magnification  of  the  trans¬ 
mitter  acousto-optic  modulator  in  the  plane  of  the  receiver  acousto- 

j 

optic  correlator.  Lg  then  provided  the  Fourier  transformation  of  j 

this  spatial  light  distribution  producing  distinct  spots  in  its 
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(a)  l*na  nodal  of  optical  waveguide 


Fig.  3.20  Optical  Coisiaunication  Syetem 
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Fig.  3.21  Optical  Conwunicationa ■  Implementation  1 
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focal  plane  corresponding  to  the  various  diffraction  orders.  A 
photodiode  collected  light  in  the  first  order  diffracted  light, 
and  electronic  processing  of  the  photodiode  signal  yielded  the 
oscilloscope  display  of  the  desired  convolution  in  the  same  manner 
as  discussed  previously. 

An  objection  to  this  arrangement  is  that  the  optical 
transmission  line  must  carry  not  only  the  modulated  light  containing 
the  information  but  the  unmodulated  light  as  well.  An  alternative 
arrangement  is  shown  in  Fig.  3.22.  In  this  configuration  only 
the  first  order  diffracted  light  from  the  first  acousto-optic 
modulator  is  permitted  to  interact  with  the  second  acousto-optic 
modulator.  The  photodiode  output  in  this  case  provides  the 
desired  convolution  without  further  electronic  processing.  Shown 
in  Fig.  3.23  are  experimental  results  using  the  latter  arrangement. 

3.3.3  Discussion 

The  discussion  of  the  real-time  convolver  is  germane  to 
the  optical  communication  system  presented,  with  two  notable 
exceptions.  First,  the  light  is  video  detected  in  the  system  of 
Fig.  3-22.  The  consequences  of  this  are  less  noise  immunity  than 
for  heterodyne  detection  and  the  relaxing  of  the  detector  frequency 
response  requirement.  The  reduced  detector  frequency  requirements 
would  enable  the  use  of  phototransistors  rather  than  the  faster 
photodiodes,  with  the  obvious  advantage  of  large  signal  gain. 

Second,  the  system  of  Fig.  3.22  does  not  produce  the 
convolution  of  the  signals  f^(t)  and  but  rather  the  square 


Fig.  3*22  Optical  Communications :  Implementation 


of  the  convolution.  This  is  not  readily  apparent  from  the  experi¬ 
mental  result  due  to  the  fact  that  a  low  frequency  oscilloscope 
was  used  to  display  the  signal,  resulting  in  rounding  of  the 
features.  That  it  is  the  square  of  the  convolution  is  easily 
seen  by  analyzing  the  system  using  the  device  model  developed 
in  Section  3.1. 

As  only  the  first  order  diffracted  light  is  passed 

by  the  aperture  stop,  the  light  distribution  incident  on  the 

acousto-optic  correlator  is  J..[V  (v  t-z-d/2)],  which  is 

1  X  & 

incident  at  the  light  frequency  (a>  +  u>  ).  Thus  the  light  ampli- 
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tude  diffracted  by  the  correlator  is  given  by 


j(w  +  w  )t 
A(6)  =  e  a  ]  J1[V1(vat-z- 


d/2)]e 


J  V  (v  t+z+d/2)sin(cj  t+k  z) 

<-  cl  " 


dz. 

(3.1*0) 


Expanding  the  complex  exponential  in  terms  of  Bessel  functions, 
the  first  order  diffracted  light  from  the  correlator  is  given  by 
j(w  +  0)  )t  r  » 

A-j^  =  e  ]  J1[V1(vat-z-d/2)]  J1[v2(vat+z+d/2 ) ]  dz  (3.1*1) 


which  in  the  limit  of  small  and  is  the  convolution  integral. 

What  is  detected  by  the  photodetector  is  the  intensity 

.  °°c  2 

ri  =  I  J  ‘Vv1(vat  -  z  -  d/2)]  J1[V2(v&t  +  z  +  d/2)]  dz  |  ,  (3.1*2) 

—00 

which  is  incident  at  the  light  frequency  (u>  +  2co  ). 

£L 

The  convolution  of  Eq.  (3.1*1)  could  be  abstracted  by 


using  optical  heterodyning  techniques,  combining  a  coherent 
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reference  with  this  light.  If  desired  a  logarithmic  amplifier  can 
be  used  to  obtain  the  proper  convolution  from  the  detected  light, 
Eq.  (3. ^2).  However,  for  purposes  of  signal  identification,  there 
is  no  signal-to-noise  advantage  in  doing  so. 

It  is  of  interest  to  compare  the  optical  communication 
system  employing  an  acousto-optic  correlator  at  the  receiver  end 
with  the  simpler  system  in  which  uncoded  serial  data  is  trans¬ 
mitted  and  received  by  video  detection.  For  simplicity  we  assume 
that  both  systems  are  characterized  by  stochastic  noise  asso¬ 
ciated  with  transmission  through  the  channel,  with  average  power 
?c<  In  addition  we  assume  that  the  photodetector  and  associated 
circuitry  at  the  reception  end  of  both  systems  can  be  characterized 
by  a  fixed  noise  power  P^,  which  in  a  crude  approximation  is  taken 
to  be  independent  of  the  incident  light  signal.  This  is  reason¬ 
able  if  the  maximum  photodetector  shot  noise  does  not  exceed  the 
thermal  noise  of  the  load  resistor. 

Since  the  channel  noise  and  detector  noise  are  uncorre¬ 
lated,  the  signal-to-noise  ratio  at  the  detector  output  for  the 
system  transmitting  uncoded  data,  in  terms  of  the  average  signal 
power  P  ,  is 


(S/N) 


uncoded 


P  +  P„ 
c  d 


(3.»t3) 


For  the  system  employing  coded  data,  designating  by 
the  ideal  process  gain  and  by  q  the  acousto-optic  deflection 
factor  of  the  correlator,  defined  as  the  ratio  of  useful  deflected 


light  to  incident  light,  the  signal-to-noise  ratio  is 


(S/N) 


coded 


^  Vs 

n  p  +  p, 

c  a 


(3.M*) 


The  deflection  factor  is  always  equal  to  or  less  than  unity.  In 
a  Bragg  deflector,  r|  can  be  1  under  maximum  deflection  conditions, 
but  for  a  Raman-Nath  device  in  which  light  is  detected  in  the  first 

p 

diffraction  order,  .339  <  n  <  0,  since  (J, )  =  .339.  Defining 

—  —  x  max 

Y  by  Y  =  P  /P^,  the  ratio  of  the  coded  to  the  uncoded  signal-to- 
noise  ratio  is  given  by 


(S/N) 


coded 


(S/N) 


uncoded 


T)  Ni(y  +  1) 
n  Y  +  1 


(3.U5) 


For  a  noisy  channel  and  fairly  efficient  acousto-optic 
deflector,  Y  »  1  and  n  Y  >>  1» 


(S/N) 


coded 


(S/N) 


uncoded 


=  N. 


(3.1*6) 


as  expected.  The  process  gain  N^for  autocorrelation  of  an  M-bit 
pseudo-random  code  is  M,  which  can  be  made  quite  large  being 
limited  by  the  correlator  time-bandwidth  product. 

In  the  event  that  the  channel  is  noiseless,  Y=  0, 


(S/N) 


coded 


(S/N) 


=  n  N. 


uncoded 


(3.1*7) 


There  is  an  advantage  in  using  the  acousto-optic  correlator  only 
if  the  code  is  sufficiently  long  to  compensate  for  light  lost  as 


a  consequence  of  using  the  acousto-optic  device.  Moreover,  if 
the  received  light  signal  is  not  buried  in  noise  due  to  the 
channel,  correlation  techniques  are  not  appropriate. 

3.*»  Image  Scanner 

3.1+.1  Introduction 

Surface  acoustic  waves  provide  a  convenient  means  of 
scanning  a  light  intensity  distribution  for  the  purpose  of 
facsimile  transmission.  The  first  reported  experiments  employed 

/  jQ  \ 

the  photoconductive  effect  along  the  surface  of  a  semiconductor. 
Allipi  et  al. demonstrated  acousto-optic  image  scanning  with 
surface  waves  using  the  top  reflection  and  back  reflection  con¬ 
figurations.  We  present  here  an  image  scanner  for  facsimile  trans¬ 
mission  of  black  and  white  transparencies  employing  the  side  entry 
technique,  and  extend  the  capability  of  that  system  to  handling 
color  transparencies. 

The  black  and  white  image  scanner  is  shown  in  Fig.  3.2l*. 
Incident  laser  light  is  spatially  modulated  by  a  one  dimensional 
transparency  which  is  imaged  into  the  LiNbO^  crystal.  A  short 
acoustic  pulse  scans  this  image  and  provides  a  light  intensity  in 
the  first  order  diffracted  light  which  replicates  the  intensity 
distribution  as  a  function  of  time.  Scanning  of  a  two  dimensional 
transparency  is  achieved  by  mechanical  translation  of  the  trans¬ 
parency  in  the  horizontal  direction  after  each  vertical  scan  of 


the  acoustic  pulse. 
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MECHANICALLY  SCANNED 


Fig.  3.24  Black  and  White  Image  Scanner 


MECHANICALLY  SCANNED 


Fig.  3.25  Color  Image  Scanner 
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Recalling  that  in  Raman-Nath  diffraction  the  diffraction 
angle  of  the  first  order  is  proportional  to  the  light  wavelength, 

X,  it  is  a  natural  extension  of  this  system  to  provide  scanning  of 
multicolor  transparencies.  The  incident  light  is  now  composed  of 
several  colors,  such  as  red,  green,  and  blue  laser  light,  and  the 
intensity  transmission  function  is  a  color  transparency.  The 
first  order  diffracted  light  exits  at  a  different  angle  for  each 
color,  as  shown  in  Fig.  3.25.  Three  photodiodes,  one  for  each 
color,  are  then  used  to  collect  the  light.  The  three  photodiode 
signals,  after  amplification,  can  then  be  multiplexed  for  trans¬ 
mission.  Because  diffraction  is  used  rather  than  color  filters 
to  separate  color  components ,  the  system  is  capable  of  color 
separation  with  high  spectral  purity. 

An  important  factor  which  characterizes  any  image 
scanner  is  resolution.  For  an  acoustic  pulse  of  time  duration 
Tp,  the  spatial  resolution  obtained  is  given  by 

H  -  Tp  va  ,  (3.U8) 

where  v  is  the  acoustic  velocity.  The  design  of  the  transducer 
used  to  launch  the  surface  acoustic  wave  provides  a  lower  limit 
for  the  pulse  duration  Tp.  For  the  interdigital  transducers  used, 
having  5  finger  pairs,  the  transducer  bandwidth  is  20%.  Thus,  the 
minimum  pulse  duration  Tp(min)  is  given  by 

IF<1o)  ~  5/fo  ’ 

where  f  is  the  transducer  center  frequency. 
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The  amount  of  light  diffracted  is  proportional  to  the 
pulse  length,  Tp.  Therefore,  as  resolution  is  improved,  the  signal 
is  decreased,  other  factors  being  equal.  For  this  reason,  in  a 
practical  scheme  one  cannot  simply  improve  resolution  without  bound 
by  going  to  higher  acoustic  frequencies.  Ultimately,  signal-to- 
noise  considerations  will  determine  the  optimal  pulse  width,  and 
therefore,  resolution.  The  signal-to-noise  ratio  at  the  photo¬ 
detector  output  is  improved  by  increasing  the  incident  light 
intensity,  but  optical  damage  to  the  LiNbO^  places  a  limit  on 
useful  intensity. 

In  addition  to  determining  resolution,  the  distance  T  v 

P  a 

also  determines  the  spot  size  of  the  diffracted  orders,  since  this 
is  the  diffraction  limiting  aperture.  If  a  lens  of  focal  length 
F  is  used  to  form  the  Fraunhofer  diffraction  pattern  of  the  diffracted 
light,  the  spot  size,  s,  in  each  order  will  be  on  the  order  of 


s  =  — 


F  \ 


T  v 
P  a 

It  is  quite  possible  that  if  Tp  is  small  enough  the  spots  corre¬ 
sponding  to  the  first  order  diffraction  for  each  color  used  may 
overlap.  In  this  case,  a  dispersive  element  such  as  a  prism  can 
be  used  to  enhance  the  spatial  separation  of  the  different  color 
components. 

3.1*. 2  Experiment 

o  o 

In  the  experiment, green  ( 51^5  A)  and  blue  (U880  A)  light 
from  an  argon  ion  laser  was  combined  in  a  beam  splitter  with  red 
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light  (6328  A)  from  a  HeNe  laser.  The  beam  was  horizontally 

expanded,  then  truncated  and  recollimated,  so  that  the  light 

incident  on  the  color  transparency  was  of  uniform  intensity  in  the 

one  dimension  of  interest.  The  surface  acoustic  waves  were  launched 

on  YZ  LiNbO^  by  interdigital  transducers  matched  to  50  ft  at  center 

'frequency.  Several  different  LiNbO^  samples  were  used,  with  center 

frequencies  ranging  from  50  to  100  MHz.  We  take  f  =  100  MHz  for 

the  following  numerical  example.  From  Eq.  (3.1*9),  Tp(m£n)  =  50  nsec. 

Using  v  =  3488  m/sec,  Eq.  (3.48)  provides  that  this  corresponds 
& 

to  a  spatial  resolution  on  the  crystal  of  0.17  mm.  This  provides 
a  scan  of  2  x  10  spots/sec.  For  a  LiUTbO^  sample  5  cm  long  the 
number  of  resolvable  spots  per  scan  is  294.  Electronic  constraints, 
however,  required  the  use  of  200  nsec  pulses.  Therefore,  the 
resolution  obtained  was  0.68  mm,  corresponding  to  73  resolvable 
spots  and  a  scan  rate  of  5  x  10^  sec  For  the  red,  green,  and 
blue  laser  beams  the  respective  diffraction  angles  for  the  first 
order  diffracted  light  were  0.019,  0.015,  and  0.014  radians  for  an 
acoustic  frequency  of  100  MHz.  Using  a  30  cm  focal  length  lens  to 
form  the  diffraction  pattern  of  the  diffracted  light  and  an  acoustic 
pulse  of  200  nsec,  the  spot  sizes  of  the  focused  diffracted  light 
were  0.27  mm,  0.22  mm,  and  0.21  mm  for  the  red,  green,  and  blue 
respectively.  These  spots  were  spatially  separated. 

Figures  3.26-3.31  show  results  obtained.  Figure  3.26(a) 
is  the  trace  of  the  signal  obtained  with  no  modulating  transparency 
present,  except  for  a  window  aperture.  The  response  shown  is 
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Fig.  3.26  Unprocessed  Scanner  Output 

(a)Photodiode  output  with  no  input  transparency  (b)  Photodiode 
output  with  input  transparency  composed  of  .8mm  wide  black 
lines  evenly  spaced  (lp  sec/div.) 
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Fig.  3.28  Two  Dimensional  Image  Scanning 
Photodiode  output  modulates  2 -axis  of . oscilloscope  while 
transparency  is  scanned  mechanically  in  the  vertical  direction 
Letters  are  4mm  high.  Acoustic  frequency=100  MHz. 


Fig.  3*29  Results  of  Color  Image  Scanning 
Photodiode  output  for  (a)  Blue  light  (b)  Red  light  (c) 
Green  light.  Optical  input  is  a  green-red-blue  striped 
transoarency .  Time  scale=lu  sec/div.  Acoustic  frequency 
=100  :  nz. 
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Fig.  3.30  Results  for  Color  Image  Scanning 
Photodiode  output  for  (a)  Blue  light  (b)  Red  light  (c) 

Green  light.  Optical  input  is  a  green-red-blue  striped 
transparency  in  which  only  one  section  is  illuminated  at 
a  time.  Time  scale-  lp  sec/div.  Acoustic  frequency  =100  MHz. 
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Fig*  3*31  Image  Scanner  Intensity  Response  (Photodide 
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output  across  5°  ohm  load) 
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k.5  psec  long  which  corresponds  to  1.6  cm  of  interaction  length. 
The  acoustic  pulse  was  200  nsec  long  at  a  center  frequency  of 
100  MHz.  Figure  3.26(b)  is  the  oscilloscope  trace  for  a  mono¬ 
chromatic  object  composed  of  .8k  mm  black  and  white  stripes. 

The  static  variations,  resulting  from  beam  nonuniformities  and 
crystal  surface  defects,  were  eliminated  by  imposing  a  threshold 
detector  and  limiter  before  display.  Results  are  shown  in  Fig. 
3.27.  The  cost  of  this  modification  is  the  digital  nature  of 
the  output  in  which  some  information  concerning  gray  tones  is 
lost . 

The  two  dimensional  scanning  capability  is  shown  in 
Fig.  3.28,  in  which  mechanical  scanning  of  the  black  and  white 
transparency  was  used  for  the  vertical  direction.  The  photodiode 
output  was  used  to  modulate  the  oscilloscope  z-axis  to  obtain 
this  picture.  A  voltage  ramp  synchronized  with  the  vertical 
mechanical  scan  was  applied  to  the  y-axis,  and  sweep  triggering 
was  synchronized  to  the  acoustic  pulse. 

Figure  3-29  demonstrates  color  scanning  operation. 

The  one  dimensional  object  was  composed  of  parallel  colored 
stripes,  green,  red,  and  blue.  The  top  trace  is  the  detector 
output  placed  to  collect  the  blue  light,  the  middle  trace  is 
for  the  red  light,  and  bottom  trace  is  for  the  green  light. 

Note  that  a  blue  component  of  light  is  in  evidence  in  the  light 
passed  by  the  green  stripe,  and  visa  versa.  Covering  the  blue 
stripe  while  detecting  the  green  light,  and  separately  covering 


the  green  stripe  while  detecting  the  blue  light,  gave  the  corre¬ 
sponding  output  signals  of  Fig.  3.30. 

Using  neutral  density  transparencies  the  photodiode  out¬ 
put  of  the  image  scanner  was  measured  as  a  function  of  light  inten¬ 
sity  incident  on  the  LiHbO^  crystal.  Results  are  shown  in  Fig.  7- 

3.1*. 3  Discussion 

The  acousto-optic  interaction  provides  a  simple  means 
of  scanning  monochromatic  as  well  as  color  transparencies  for 
the  purpose  of  transforming  the  light  intensity  distribution 
into  a  voltage  signal.  In  the  system  demonstrated,  resolution 
was  limited  by  the  bandwidth  of  a  preamplifier  in  the  detection 
system,  but  normally  resolution  would  be  limited  by  transducer 
bandwidth.  Although  higher  bandwidths  can  be  achieved  with 
bulk  wave  transducers,  the  higher  acoustic  velocity  of  bulk 
waves  tends  to  degrade  resolution.  For  this  reason  it  is 
advantageous  to  use  surface  acoustic  wave  devices  for  image 
scanning. 

The  lens  system  used  in  imaging  the  transparency  into 
the  crystal  can  be  used  to  provide  magnification.  Hence  resolu¬ 
tion  is  a  measure  of  spatial  separation  achieveable  in  the  plane 
of  the  crystal,  not  in  the  transparency  plane.  For  this  reason 
a  better  measure  of  system  capability  is  the  number  of  resolveable 
spots.  The  number  of  resolveable  spots  is  proportional  to  the 
transducer  bandwidth  and  acoustic  propagation  length,  hence  it  is 


directly  proportional  to  the  device  time-bandwidth  product.  It  is 

profitable  to  improve  the  time-bandwidth  product  by  means  of 

increasing  transducer  bandwidth,  rather  than  going  to  longer 

crystals,  because  for  two  dimensional  scanning  the  minimum  time 

for  each  individual  line  scan  is  given  by  the  transit  time  of  the 

acoustic  pulse  to  travel  the  sample  length. 

Because  laser  sources  were  used  to  illuminate  the 

transparency,  the  color  components  detected  are  not  directly 

compatible  with  commercial  television  standards.  The  HeNe  red 

is  only  slightly  deeper  than  that  specified  by  the  FCC.  The 

blue  and  green  argon  lines,  however,  depart  considerably  from 

both  the  FCC  values  and  the  currently  used  television  phosphors, 
o  o  o 

The  51^5  A  green  should  be  closer  to  5350  A,  and  the  U880  A  blue 

o 

should  be  1+700  A.  The  use  of  the  argon  lines  excludes  discrimina¬ 
tion  of  deep  blues  and  purples.  It  is  to  be  noted  that  the 

imaging  system  is  completely  compatible  with  the  laser  TV  display 

(79) 

system  reported  by  Korpel  et  al .  ,  as  they  employ  the  same 
laser  lines.  No  colorimetric  distortion  would  be  introduced  in 
using  the  two  complementary  systems,  but  as  pointed  out  deep 
blues  and  purples  would  be  absent. 

An  especially  attractive  feature  of  the  acousto-optic 
image  scanner  is  that  it  is  easily  adapted  to  other  functions 
such  as  electronic  focusing  of  images  and  forming  real-time  Fourier 
transforms  of  images.  These  operations  are  achieved  electronically 
without  any  alteration  of  the  optical  system.  For  example,  to 
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obtain  the  Fourier  transform  of  a  transparency  imaged  into  the 
LiNbO^,  two  acoustic  waves  are  launched  from  opposite  ends  of 
the  crystal.  One  is  the  usual  narrow  scanning  pulse  which  modu¬ 
lates  its  carrier  which  is  at  frequency  to  ,  while  the  second  is  a 

Q. 

continuous  constant  amplitude  wave  at  the  frequency  to  +  Ato  . 

EL  & 

The  photodiode  voltage,  processed  as  for  the  ambiguity  function, 
provides  the  signal  t(v  t)  exp  [Ato  t].  Integration  of  this 

cL  8. 

signal  in  the  time  domain  provides  the  desired  Fourier  transform 
as  Ato^  is  swept  in  time. 


CHAPTER  k 


SUMMARY  AND  CONCLUSIONS 

l+.l  Summary 

A  classical  description  of  the  interaction  of  laser  light 
with  ultrasonic  surface  acoustic  waves  and  an  experimental  demon¬ 
stration  of  the  application  of  the  interaction  to  several  new 
devices  is  presented.  The  analysis  assumes  linear  constitutive 
relations,  and  it  assumes  the  perturbation  of  the  impermeability 
tensor  elements  is  linearly  related  to  the  strain  and  electric 
field  associated  with  the  acoustic  wave.  The  inclusion  of  a  term 
proportional  to  the  rotation  tensor  is  readily  incorporated  into 
the  analysis  but  is  deleted  in  the  numerical  calculation  due  to 
the  unavailability  of  published  data  for  the  appropriate  coeffi¬ 
cients.  Good  agreement  with  previously  published  experimental 
data  is  achieved  but  requires  some  modification  of  the  reported 
photoelastic  constants. 

A  scalar  diffraction  model  of  the  acousto-optic  inter¬ 
action  is  discussed  and  is  used  in  describing  the  operation  of  the 
modulator,  image  scanner,  and  convolver.  The  scalar  diffraction 
model  is  related  analytically  to  the  theoretical  treatment  of  the 
interaction  in  surface  acoustic  waves  and  numerically  to  experi¬ 
mental  results  for  the  convolver. 

New  experimental  results  are  obtained  and  device  appli¬ 
cations  demonstrated  using  the  side-entry  interaction  configuration. 
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New  experimental  results  for  the  real-time  convolver  employing 
optical  heterodyning  include  autocorrelation  of  Barker  codes, 
ambiguity  function  generation,  and  the  measurement  of  diffraction 
characteristics  for  y-  and  z-polarized  light. 

Results  are  given  for  the  correlation  of  a  time  varying 
light  intensity  distribution  with  an  electrical  signal,  and  the 
application  to  an  optical,  communication  system  is  discussed.  A 
scheme  involving  an  acousto-optic  modulator  at  the  transmission 
end  and  an  acousto-optic  correlator  at  the  reception  end  is 
demonstrated.  In  this  scheme  only  modulated  light  in  the  first 
diffraction  order  of  the  acousto-optic  modulator  is  transmitted 
with  the  result  that  squared  correlation  is  performed  at  the 
reception  end,  detected  in  real-time  by  video  detection. 

The  scanning  of  one  and  two  dimensional  black  and  white 
transparencies  is  performed  using  the  side-entry  technique. 
Mechanical  scanning  is  employed  for  the  two  dimensional  scan,  but 
could  be  replaced  by  faster  light  deflection  methods.  The  system 
is  extended  to  the  scanning  of  color  transparencies.  High  spectral 
purity  is  achieved,  and  no  cross-talk  due  to  overlapping  diffrac¬ 
tion  orders  is  observed. 

h.2  Conclusions 

The  side-entry  technique  provides  a  simple  and  efficient 
means  of  diffracting  an  unguided  optical  beam  by  means  of  the 
acousto-optic  interaction  using  surface  acoustic  waves.  The 
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technique  provides  a  much  longer  interaction  length  than  any  of 
the  other  interaction  configurations,  comparable  to  bulk  wave 
devices,  while  retaining  the  advantageous  features  of  surface 
acoustic  waves .  These  features  include  high  power  density  for 
a  given  acoustic  power,  a  lower  acoustic  group  velocity,  and 
use  of  planar  fabrication  techniques.  In  using  Rayleigh  waves, 
there  is  the  additional  significant  advantage  that  the  acoustic 
propagation  is  dispersionless.  The  potential  of  utilizing 
thin  films  is  also  available  with  surface  wave  devices  which  is 
not  with  bulk  wave  devices.  The  longer  interaction  length 
achieved  with  the  side-entry  technique  which  is  not  possible  in 
the  other  interaction  configurations  provides  a  more  efficient 
diffraction,  in  terms  of  electrical  drive  power,  and  puts  the 
interaction  in  the  Bragg  regime  at  a  lower  frequency.  At  1+5  MHz 
about  50%  of  the  zero  order  light  can  be  modulated  with  300  mW 
of  applied  electrical  power  using  the  side-entry  technique.  Any 
of  the  other  three  configurations  would  provide  only  a  few 
percent  modulation.  The  side-entry  configuration  can  provide 
Bragg  scattering  at  frequencies  below  150  MHz.  The  top  trans¬ 
mission  configuration  remains  in  the  Raman-Nath  regime  to  fre¬ 
quencies  above  1  GHz. 

In  terms  of  electrical  drive  power,  the  devices  reported 

here  employing  the  side-entry  technique  do  not  compare  favorably 

with  the  best  integrated  optics  acousto-optic  modulators  reported 

( 57 ) 

in  recent  literature.  Nguyen  reports  acousto-optic  Bragg 


deflection  employing  guided  optical  waves  and  surface  acoustic 
waves  in  which  50$  of  the  output  zero  order  light  is  deflected 
with  68  mW  of  electric  drive  power.  This  result  can 
probably  only  be  achieved  with  integrated  optics  structures. 

The  paper  cited,  however,  does  not  list  optical  losses  associated 
with  coupling  into  the  waveguide,  so  actual  fractional  deflection 
might  be  less  than  50$. 

All  of  the  devices  experimentally  reported  here  would 
benefit  by  increasing  the  time-bandwidth  product.  To  use  spread 
spectrum  techniques  to  advantage  in  a  communications  system 
employing  real-time  correlation  at  the  receiver  it  is  necessary 
to  handle  codes  of  long  bit  length  to  obtain  substantial  process 
gain.  At  best,  the  process  gain  using  signal  correlation  to 
identify  a  signal  comprised  of  M  bits  is  M,  except  for  the  case 
of  Barker  codes  in  which  M2  is  achieveable  for  M  £  13.  Because 
the  operation  of  the  acousto-optic  convolver  requires  that  all 
M  bits  must  simultaneously  be  accommodated  on  the  device  acoustic 
propagation  path  for  maximum  correlation  to  occur,  the  crystal 
sample  must  be  long  or  the  spatial  extent  of  each  bit  must  be 
short.  Moreover,  the  maximum  code  length  is  proportional 
to  the  device  time-bandwidth  product.  In  the  case  of  ambiguity 
function  generation  for  radar  or  sonar,  discrimination  of  a 
target's  range  and  velocity  is  proportional  to  the  time-bandwidth 
product.  In  image  scanning  applications  of  the  acousto-optic 
interaction,  the  number  of  resolveable  spots  per  scan  is  limited 
by  the  device  time-bandwidth  product. 
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It  is  readily  apparent  that  for  the  applications  reported 
the  transducer  bandwidth  should  be  large,  and  the  optical  aperture 
should  be  wide,  except  as  noted  for  the  image  scanner  if  a  very 
high  repetition  rate  is  required  for  two  dimensional  scanning. 
Generally  it  is  impractical  to  use  an  optical  aperture  greater 
than  ^  6  cm,  although  crystals  of  longer  dimension  are  readily 
available.  It  might  be  possible  to  employ  a  folded  acoustic  path 
to  increase  the  time  bandwidth  product  of  the  convolver,  but  this 
method  would  not  be  of  use  for  the  image  scanner.  The  only  other 
means  of  obtaining  higher  time-bandwidth  products  is  to  use  wider 
transducer  bandwidths.  For  fixed  fractional  bandwidth  this 
requires  higher  frequency  operation.  Effective  bandwidth  can  be 
increased,  also,  by  using  a  multiple  transducer  array  composed 
of  transducers  of  sequential  center  frequencies  with  overlapping 
or  contiguous  bandwidths.  With  present  technology  bandwidths 
of  several  hundred  megahertz  are  attainable,  providing  potential 
time-bandwidth  products  of  several  thousand.  A  6  cm  optical 
aperture  and  a  200  MHz  bandwidth  would  provide  a  time-bandwidth 
product  of  3^0,  which  is  well  within  the  capability  of  present 
transducer  technology.  It  is  significant  that  such  time-bandwidth 
products  are  not  to  be  expected  from  surface  acoustic  wave  inte¬ 
grated  optics  devices  because  of  difficulties  encountered  in 
coupling  over  large  optical  apertures.  A  1  cm  aperture  is  the 
largest  reported  to  date.^®^  Of  course  improvement  is  to  be 


expected. 
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the  first  order  intensity  for  Bragg  diffraction  to  that  for  Raman- 

1  2 

Nath  diffraction,  [sin  x)/J^(x)]  ,  as  a  function  of  the  argu¬ 
ment  x.  Thus  for  most  applications,  there  is  no  disadvantage  in 
going  to  the  higher  acoustic  frequencies  at  which  the  Bragg 
regime  occurs. 

2 

sin(ix) 

.  «v*>. 
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Fig.  4.1  Comparisoh  of  Bragg  and  Raman-Nath  First  Order 

Intensity 


There  is  one  major  difficulty  attendant  with  Bragg 
diffraction,  and  that  is  the  optical  aperture  characteristic. 
For  a  given  light  wavelength  and  angle  of  incidence,  efficient 
diffraction  occurs  only  for  a  very  small  bandwidth. 
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It  is  of  interest  that  while  wider  bandwidths  are  avail¬ 
able  using  bulk  acoustic  waves,  this  does  not  result  in  a  propor¬ 
tionate  gain  in  time-bandwidth  product.  The  greater  bulk  wave 
velocity  produces  a  diminution  of  time-bandwidth  product. 

As  was  indicated  in  Chapter  1,  the  acoustic  aperture 
of  an  interdigital  transducer  is  important  in  determining  input 
impedance.  For  a  four  and  one  half  finger  pair  transducer,  which 
optimizes  the  efficiency-bandwidth  product  for  a  transducer 
fabricated  on  YZ-LiNbO^,  the  acoustic  aperture  which  yields  a 
50  ft  resistive  component  of  input  impedance  is  108  A.  Since  the 
interaction  length  in  the  side-entry  configuration  is  given  by 
the  acoustic  aperture,  this  has  the  interesting  consequence  that 
if  impedance  matching  is  to  be  achieved  without  the  use  of  a 
transforming  network,  the  transducer  center  frequency  determines 
whether  the  interaction  is  to  be  in  the  Raman-Nath  or  Bragg 
regime.  For  frequencies  of  100  MHz  and  below,  the  interaction 
is  in  the  Raman-Nath  regime.  For  frequencies  in  excess  of 
100  MHz  the  interaction  becomes  more  Bragg-like.  This  is  an 
important  consideration. 

For  low  acoustic  powers,  the  Raman-Nath  and  Bragg 
regimes  yield  the  same  amount  of  light  deflected  into  the  first 
diffraction  order.  This  follows  from  Eqns.  (2.65)  and  (2.77), 
using  the  approximations  sin  x  ^  x  and  J^(2x)  'v  x.  However  at 
higher  powers,  the  Bragg  interaction  is  more  efficient.  This  is 
shown  by  the  graph  below.  Fig.  4.1,  which  shows  the  ratio  of 


frequency  which  satisfies  the  Bragg  condition.  No  such  difficulty 
is  associated  with  Raman-Nath  devices.  Hence,  signal  processors 
in  the  Bragg  regime  require  phased-array  frequency-steering  acoustic 
transducers  for  wide  band  operation. 

In  the  color  image  scanning  experiment,  three  color 
o  o 

components  ranging  from  1*880  A  to  6328  A  are  incident  in  a  colli¬ 
mated  beam  on  the  LiNbO^  crystal.  Only  Raman-Nath  diffraction 
provides  a  simple  means  of  deflecting  all  the  colors  with  the 
same  acoustic  pulse;  Bragg  diffraction  is  inappropriate  for  this 
application.  As  a  consequence,  at  very  high  acoustic  frequencies 
necessary  to  achieve  a  large  number  of  resolveable  spots,  the 
side-entry  technique  is  not  suitable.  One  of  the  less  efficient 
interaction  configurations,  such  as  top-transmission,  must  be 
used  in  order  to  remain  in  the  Raman-Nath  regime. 

It  should  be  pointed  out  that  as  the  acoustic  fre¬ 
quency  increases,  the  acoustic  power  resides  proportionally 
closer  to  the  crystal  surface.  One  can  anticipate  that  guiding 
of  the  optical  wave  to  confine  it  near  the  crystal  surface  will 
be  essential  at  high  acoustic  frequencies.  For  this  reason  it 
is  to  be  expected  that  most  research  in  the  near  future  in  surface 
wave  acousto-optics  will  be  in  the  area  of  integrated  optics. 

From  the  foregoing  discussion,  the  following  guidelines 
concerning  the  application  of  the  acousto-optic  interaction  in 
surface  acoustic  waves  seem  apropos.  From  the  point  of  view  of 
power  consumption,  integrated  optics  devices  are  most  attractive 


and  will  be  even  more  so  as  higher  acoustic  frequencies  are 
achieved.  These  devices  are  especially  well  suited  at  the 
present  time  for  such  functions  as  digital  switching,  modula¬ 
tion,  and  beam  steering,  since  a  narrow  optical  beam  can  be 
used  in  these  operations.  For  device  applications  in  which  a 
large  time  bandwidth  product  is  required,  such  as  those  appli¬ 
cations  discussed  in  this  report,  the  wider  optical  aperture 
available  with  the  side-entry  technique  can  be  used  to  advan¬ 
tage.  Should  large  time-bandwidth  products  be  obtained  with 
integrated  optics  devices  in  the  future,  almost  all  the  device 
applications  demonstrated  here  using  the  side-entry  technique 
could  be  performed  equally  well  using  a  guided  optical  wave. 

The  one  exception  is  that  of  the  color  image  scanner.  In  order 
to  remain  in  the  Raman-Nath  regime  at  higher  acoustic  frequen¬ 
cies,  it  will  be  necessary  to  employ  a  less  efficient  inter¬ 
action  configuration. 

It  is  of  interest  to  point  out  several  specific  areas 
where  further  investigation  of  the  acousto-optic  interaction  in 
surface  acoustic  waves  would  be  pertinent.  An  analysis  of  the 
effect  of  cone  angle  and  incident  light  inclination  is  espe¬ 
cially  relevant  to  the  side-entry  technique.  In  the  present 
treatment  a  truncated  uniform  plane  wave  with  propagation  in  the 
x-z  plane  is  assumed,  but  this  idealization  does  not  obtain  in 


practice. 


An  experimental  program  to  determine  the  coefficients  of 


the  rotation  contribution,  would  be  of  significant  value. 

Such  coefficients  might  account  for  the  discrepancies  in  the 
literature  for  the  other  photoelastic  constants. 

In  the  area  of  applications,  the  present  work  suggests 
a  number  of  new  applications  and  new  devices  which  should  be 
investigated.  In  particular,  in  the  signal  processor  described 
in  Chapter  3,  the  integration  associated  with  computing  the  con¬ 
volution  is  performed  in  the  spatial  domain.  In  many  applications 
in  which  careful  laboratory  optical  alignment  is  not  practical, 
it  might  be  desireable  to  perform  this  integration  in  the  time 
domain.  Implementation  involves  using  a  narrow  optical  aperture 
and  adding  a  low  pass  filter  for  the  final  electronic  processing 
stage. 


The  same  idea  is  readily  carried  over  to  other  device 
and  system  applications.  The  optical  communication  system  dis¬ 
cussed  in  the  previous  chapter  requires  that  the  transmission 
channel  preserve  the  amplitude  distribution  of  the  optical  wave- 
front.  While  this  may  not  be  a  deficiency  for  an  intersatellite 
communication  network,  it  is  far  too  stringent  a  requirement  for 
present  day  optical  waveguides  used  in  conventional  terrestrial 
communications.  An  alternative  system  which  does  not  impose  this 
requirement  could  again  be  implemented  by  taking  a  narrow  optical 
beam,  which  we  may  consider  a  spatial  <5- function,  and  using  a  low 
pass  filter  to  facilitate  temporal  integration.  In  such  a  system 


the  transmitter  end  need  not  be  an  acousto-optic  device.  Any 
suitable  means  of  light  modulation  could  be  used  for  transmission, 
including  direct  modulation  of  a  solid  state  emitter.  The  impor¬ 
tant  feature  of  such  a  system  is  the  acousto-optic  correlator  at 
the  receiver  end. 

A  slight  variant  of  the  optical  communication  system 
described  above  could  be  used  for  pulse-doppler  laser  radar.  In 
this  case  a  high  power  laser  and  bulk  wave  acousto-optic  modulator 
would  be  required  at  the  transmission  end.  Following  acousto¬ 
optic  correlation  at  the  receiver,  heterodyne  optical  detection 
to  determine  doppler  shift  is  easily  performed  since  transmitter 
and  receiver  can  be  located  in  close  proximity,  thus  allowing 
the  use  of  the  same  laser  for  purposes  of  local  oscillation. 

Another  area  which  warrants  attention  is  the  use  of 
the  collinear  interaction  for  signal  processing  purposes.  Because 
polarization  rotation  occurs  only  over  a  narrow  acoustic  bandwidth 
for  a  given  light  wavelength,  previous  attention  has  been  paid 
only  to  narrow  band  filtering.  Barker  encoded  data,  however, 
imposed  upon  a  carrier's  phase,  can  be  accommodated  with  narrow 
band  signals.  The  optical  arrangement  for  a  correlator  in  such 
a  system  is  shown  in  Fig.  k.2.  The  desired  correlation  signal 
appears  in  the  photodiode  output  at  twice  the  acoustic  frequency. 
By  using  a  GaAs  LED  for  the  light  source  and  a  compatible  silicon 
p-i-n  photodiode  for  detection,  operation  using  YZ-LiNbO^  would 
be  in  the  300  MHz  range.  No  lenses  are  required  for  optical 


alignment,  and  the  crystal  itself  would  provide  directional  guiding 
of  the  LED  output.  Simplicity  of  alignment  and  compact  packaging 
contribute  to  the  commercial  attractiveness  of  the  device. 

As  with  all  devices  employing  the  side-entry  configura¬ 
tion,  there  is  a  corresponding  integrated  optics  analog.  An 
integrated  optics  collinear  convolver  could  be  built  based  on 
waveguide  mode  conversion,  rather  than  polarization  rotation. 

Again,  operation  is  narrow  band,  but  the  frequency  of  operation 
is  determined  by  waveguide  thickness,  which  is  controlled  by  the 
designer.  In  the  device  utilizing  polarization  rotation,  the 
operating  frequency  is  dictated  by  the  crystal's  birefringence. 

The  achieveable  performance  of  devices  employing  the 
acousto-optic  interaction  in  surface  waves  is  intimately  tied 
to  surface  acoustic  wave  technology  in  general.  Progress 
including  improvements  in  transducer  technology  and  the  identi¬ 
fication  of  new  crystal  cuts  with  desireable  properties  will 
undoubtedly  be  reflected  in  more  efficient  acousto-optic  devices 
operating  at  higher  frequencies  and  wider  bandwidths.  Further 
improvements  will  come  as  a  result  of  progress  in  acousto¬ 
optic  materials  research.  It  is  reasonable  to  expect  that 
surface  wave  acousto-optic  technology,  judiciously  applied, 
will  result  in  economical  devices  to  augment  the  existing  options 
of  the  optical  and  system  engineer. 
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DERIVATION  OF  THE  WAVE  EQUATION 

In  this  appendix  we  derive  from  Maxwell's  equations 
of  electromagnetic  theory  the  wave  equation  which  is  the 
starting  point  for  the  analysis  of  the  acousto-optic  inter¬ 
action.  In  rationized  MKS  units.  Maxwell's  equations  can  be 
written  in  the  form 


V  •  e  E 

=  P 

(A.l) 

V  x  E  = 

3B 

at 

(A. 2) 

V  •  B  = 

0 

(A. 3) 

-»  a  * 

(A.h) 

V  x  B  = 

y  J  +  -XT-  e  E 

Ko  at 

where  it  has  been  assumed  that  the  material  is  linear  and  non- 

magnetic.  E  and  B  are,  respectively,  the  electric  field  and 

* 

magnetic  field  intensity,  p  and  J  are  the  free  charge  density 

and  free  current  density,  is  the  permeability  of  free  space, 

** 

and  e  is  the  permittivity  tensor  which  describes  the  material. 

The  equation  of  interest  is  obtained  from  the  vector  equation 
(A. 2)  in  the  usual  manner  by  taking  the  curl  of  that  equation, 

A  2  ^  4 

using  the  vector  identity  VxVxA  =  -V  A  +  V(V  •  A)  for 
any  vector  A,  reversing  the  order  of  the  curl  and  time  derivative 
operations  on  the  right  hand  side,  eliminating  B  by  means  of 
Eq.  (A.4),  and  assuming  p  =  J  =  0  in  the  material. 
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Expressing  the  resultant  vector  equation  in  an  ortho- 


1U5 


1 


gonal  coordinate  system  in  which  E  *  e^  E^,  where  e^  is  a  unit 
vector,  and  the  Einstein  convention  is  employed  in  which  repeated 
subscripted  indices  indicate  summation  over  those  indices,  we 
obtain  the  three  scalar  equations: 


V2  E.  +  [7(7  •  E)].  =  -  u-$—  e.  E, 
i  0  3t2  ^  J 


(A. 5) 


Summation  over  j  =  1,2,3  is  understood  on  the  right  hand  side; 
the  subscript  i  can  take  on  values  i  =  1,2,  or  3.  We  shall  now 
show  that  for  all  cases  which  we  shall  consider,  the  term  on  the 
left  hand  side  involving  the  gradient  is  negligible  compared  to 
all  other  terms  in  the  equation.  We  indicate  this  for  the 
isotropic  case  in  which  the  permittivity  can  be  represented  by 
a  scalar  function,  however  the  result  will  apply  to  a  tensor  as 
well. 

With  p  =  0,  Eq.  (A.l)  can  be  written 


7  * 


E 


(A. 6) 


in  obtaining  which  the  following  vector  identity  has  been 
employed: 


7 


f  F 


F  •  7  f  +  f  7  •  F 


(A. 7) 


in  which  f  is  a  scalar  function  and  F  is  a  vector-valued  function 
of  a  scalar  argument.  In  Eq.  (A. 6),  E  is  taken  to  be  the  electric 
field  associated  with  light  propagation,  and  e  represents  the 
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permittivity  of  the  material,  which  we  allow  to  be  perturbed  by 
the  presence  of  a  propagating  acoustic  wave.  For  light  propagating 
in  the  x^-direction,  or  nearly  so,  we  may  neglect  and  take 

E  E2  +  E3  (x,y,z,t)  ,  (A. 8) 

in  which  the  orthogonal  axes  {x^,  x x^)  are  identified  with 
the  usual  (x,y,z)  designations.  Thus  the  inner  product  on  the 
right  hand  side  of  Eq.  (A. 6)  only  involves  a  sum  over  the  y 
and  z  components. 

It  will  be  sufficient  to  consider  an  exponentially 
decaying  surface  acoustic  wave  propagating  in  the  z-direction. 

In  this  case  the  relative  permittivity  can  be  written  as 

~  e  =  e°  +  Ae”ay  cos  (a>  t  -  k  z)  (A. 9) 

t  «•  » 

o 

where  eQ  is  the  permittivity  of  free  space,  e°  is  the  unperturbed 
relative  permittivity,  and  the  subscript  "a"  refers  to  the 
acoustic  wave.  As  the  perturbation  is  very  small,  A/e°  <  10~  , 
we  may  take  e  ^  Cq  e°  in  the  denominator  of  Eq.  (A. 6),  and  write 
(A. 5)  in  the  isotropic  case  as 


-,2  Ei  -  a!:  (E  '  *§>  - 
i  e 


— \  -^-r  {[e°  +  A  e  ay  cos  (co  t  -  k  z)]  E. } 
c2  9t2  a  a  i 


(A. 10) 


where  we  have  used  the  fact  that  the  speed  of  light  in  free  space, 
c,  is  given  by  (eq  pQ)  .  In  the  anisotropic  case,  in  which  we 
allow  e  to  be  a  tensor,  the  last  term  of  the  right  hand  side  would 
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be  replaced  by  the  sum  of  two  terms ,  each  of  which  would  be 

similar  in  form  to  the  term  given  in  Eq.  (A.IO). 

We  show  now  that  the  gradient  term  on  the  left  hand 

side  of  Eq.  (A.IO)  is  negligible  compared  to  the  other  small 

term  in  the  equation,  -l/c^  9^/3t^  A  e  cos  (u)  t  -  k  z)  E. . 

&  &  1 

In  establishing  this  result,  we  make  the  physically  reasonable 
assumption  that  the  spatial  rate  of  change  of  the  electric 
field  E^,  normalized  to  the  magnitude  of  E^ ,  is  less  than  or 
equal  to  that  of  the  permittivity  £. 

The  relation  we  wish  to  establish  is 


Z 

J-2,3 


A  e_Cty  cos  (id  t  -  k  z)  E.  | 

8L  &  1 


«  l 


(A. 11) 


Consider  one  term  of  the  numerator  of  (A. 11).  Using  the  assump¬ 
tion  that  the  spatial  rate  of  change  of  the  electric  field  is 
less  than  that  of  the  permittivity,  observing  that  the  surface 
wave  decay  constant  a  is  typically  on  the  order  of  k  ,  and 
making  use  of  Eq.  (A. 9), 


-J_  (JL 
9XJ  V 


)  I 


< 
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A  e 


-ay 


cos 


(ti)  t 
a 


-  k  z) 
a 


(A. 12) 
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The  optical  electric  field  is  periodic  in  time  with 

an  angular  frequency  u  which  is  much  larger  than  to  .  Thus  the 

& 

denominator  of  (A. 11)  obeys 


t  a2 

I  ~2  — 2  e'<Xy  cos  (“a*  "  ka2^  Ei^ 

c  3t 

>  |  a)2  A  e"07  cos  (wat  -  k&z)  Ei  |  (A. 13) 

c 

The  relation  to2/ c2  »  k  2,  together  with  Eqs.  (A. 12)  and  (A. 13), 
establishes  the  inequality  (A.ll). 

It  is  apparent  that  the  same  arguments  could  easily 
be  generalized  to  the  anisotropic  case.  The  final  result  is 
that  the  wave  equation  (A. 5)  can  be  written  neglecting  the 
gradient  term: 


(A.  lit) 


